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Chapter  1 
INTRODUCTION 

Although  the  elastic  finite-displacement  response  of  shallow  cylin¬ 
drical  ly  curved  panels  has  been  subjected  to  extensive  consideration  in 
the  past  fifty  years,  very  little  work  has  been  done  on  the  behavior  of 
such  panels  fabricated  from  orthotropic  materials.  The  utilization  of 
composite  materials  in  primary  structural  components  in  aircraft  appli¬ 
cations  is  steadily  increasing.  The  current  trend  is  the  implementation 
of  composite  materials  into  the  fuselage  structure.  Figure  1  shows  how 
the  fuselage  skin  may  be  represented  by  a  shallow  cylindrically  curved 
panel . 

One  concern  that  must  be  addressed  involves  the  crashworthiness  of 
a  composite  fuselage.  Crashworthiness  studies  imply  high  speed,  large 
deflection  behavior  of  structural  components.  It'is  the  intent  of  this 
study  to  lay  the  groundwork  for  understanding  the  dynamic,  elastic, 
finite-displacement  behavior  of  shallow  cylindrically  curved  panels 
fabricated  from  a  fibrous  composite  material.  This  will  be  accomplished 
by  presenting  a  complete  examination,  both  theoretical  and  experimental, 

of  the  static  behavior  of  such  a  structure  subjected  to  a  line  load 

/' 

acting  along  a  cylinder  generator.  The  line  load  acts  toward  the  center 
of  curvature  of  the  panel.  This  examination  of  static  behavior  is  a 
necessary  step  before  any  attempts  at  determining  dynamic  behavior  are 


made 


Previous  Work 

In  1929,  Biezeno  [1]  investigated  the  response  of  a  shallow  circu¬ 
lar  arch  with  pinned  ends  at  a  fixed  distance  apart  subjected  to  a 
concentrated  load  acting  at  midspan.  Using  two  solution  methods,  one  a 
Shallow  curved-beam  solution  and  the  other  a  series  solution,  Biezeno 
presented  a  detailed  analysis  of  the  syntnetric  response.  Herein,  sytitne- 
:  try  is  defined  with  respect  to  displacement  and  the  midspan  location  on 
the  curved  member.  Figure  2a  illustrates  what  is  referred  to  as  a 
symmetric  response  while  Fig.  2b  illustrates  an  asymmetric  response.  No 
results  were  presented  by  Biezeno  concerning  the  asymmetric  response 
except  to  mention  its  existence  and  its  importance  to  loss  of  stability 
of  the  symmetric  equilibrium  configuration.  Following  the  approach  of 
Biezeno,  Fung  and  Kaplan  [2]  conducted  a  complete  investigation  into  the 
behavior  of  shallow  arches,  specifically  the  case  of  an  initially  sinu¬ 
soidal  arch,  using  a  series  solution.  They  investigated  the  asymmetric 
displacement  mode  in  greater  detail  than  did  Biezeno.  It  was  found  that 
on  the  primary  load-deflection  equilibrium  path,  corresponding  to  sym¬ 
metric  deformations,  asymmetric  equilibrium  states  occurred  at  unstable 
bifurcation  points.  Schreyer  and  Masur  [3]  performed  a  similar  analysis 
on  the  circular  arch  with  clamped  ends.  In  addition  to  determining  the 
equilibrium  behavior  of  the  structure,  they  also  performed  an  in-dept^ 
analysis  of  the  stability  of  these  equilibrium  states.  In  their  study 
the  initial  rise-to-span  ratio  determined  whether  the  buckling  of  the 
arch  occurred  at  a  limit  point  on  the  primary  path  or  whether  it 
occurred  at  the  bifurcation  point  on  the  primary  path.  These  first 
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(a)  Symmetric 


FIG.  2  Symmetric  and  asymmetric  response 
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studies  into  the  finite-deformation  behavior  of  the  shallow  arch  were 
restricted  to  isotropic  and  homogeneous  materials  with  linear  elastic 
behavior  and  did  not  consider  load  eccentricity  or  geometric  imperfec¬ 
tions. 

In  recent  years  it  became  obvious  that  a  theory  which  did  not  take 
into  account  initial  imperfections  of  the  structure  would  not  be  able  to 
accurately  model  a  real  structure.  Due  to  the  imperfection  sensitivity 
of  some  arch  structures,  it  was  at  times  difficult  to  experimentally 
verify  the  analyses  put  forward.  The  buckling  loads,  as  well  as  the 
entire  character  of  the  response,  are  affected  by  the  imperfections 


present.  Schreyer  presented  an  analysis  of  the  effects  of  initial 
imperfections  on  the  response  of  shallow  circular  arches  v/ith  clamped 
ends.  In  the  study,  imperfections  in  the  geometry  of  the  arch,  which 
had  the  same  form  as  the  buckling  mode,  and  eccentricities  in  the  load 
position  were  both  examined.  Both  types  of  imperfections  were  seen  to 
cause  a  decrease  in  the  buckling  load  for  certain  ranges  in  geometry. 
This  implied  that  the  arch  is  an  imperfection-sensitive  structure. 

I 

Plai^'t  [5]  performed  a  more  thorough  investigation  into  the  effects  of 


load  position  on  the  magnitude  of  the  buckling  load.  Examining  the 
pinned-ended  shallow  circular  arch,  a  number  of  different  behaviors  were 
observed  due  to  load  offsets.  However,  only  those  arches  which  exhibi¬ 
ted  loss  of  stability  at  bifurcation  points  when  the  load  was  exactly  at 


midspan  were  seen  to  be  sensitive  to  load  offset.  In  such  an  arch  a 
small  perturbation  in  the  load  position  could  cause  significant  reduc¬ 
tions  in  the  predicted  buckling  loads.  Cheung  and  Babcock  [6]  performed 


an  experimental  investigation  into  the  effects  of  load  offset  on  the 
static  buckling  load  of  shallow  circular  arches.  Using  aluminum  arches 
with  clamped  ends,  buckling  loads  were  measured  for  load  offsets  of  up 
to  6%  of  the  arc  length.  The  results  of  the  experiment  show  significant 
reductions  in  the  magnitude  of  the  buckling  load  with  increasing  load 
offset.  Taking  into  account  such  load-position  effects  considerably 
improved  the  understanding  of  the  response  of  the  shallow  arch.  How¬ 
ever,  the  vast  majority  of  the  work  has  concentrated  on  metallic 
arches.  With  composite  materials  playing  an  increasingly  important  role 
in  the  aerospace  industry,  the  effects  of  orthotropic  materials  on  the 
response  of  the  shallow  arch  was  still  needed. 

Some  work  on  the  behavior  of  orthotropic  shallow  arches  has  been 
done.  Nash  and  Hsu  [7]  analyzed  the  structural  behavior  of  a  shallow 
circular  arch  fabricated  from  a  composite  n^aterial  with  the  fiber  direc¬ 
tion  normal  to  the  middle  surface.  In  their  theory,  Nash  and  Hsu  inclu¬ 
ded  the  effects  of  transverse  shear  stresses  and  transverse  normal 
stresses.  The  qualitative  behavior  of  their  orthotropic  arch  was  essen¬ 
tially  the  same  as  for  an  isotropic  arch.  As  is  the  case  with  an  iso¬ 
tropic  arch,  they  found  the  behavior  of  the  arch  was  controlled  by  a 
geometric  parameter.  In  addition  to  the  analytic  considerations,  they 
also  conducted  experiments  into  the  load-displacement  relationship  up  to 
the  buckling  load.  In  the  results  obtained,  a  fair  agreement  between 
theory  and  experiment  was  seen.  Marshall,  Rhodes,  and  Banks  [8]  conduc¬ 
ted  experiments  into  the  behavior  of  shallow  curved  composite  panels, 
also  finding  good  agreement  between  theory  and  experiment.  These  exper- 


f 


.  V  '  / 


7 

iments  were  performed  on  spherically  curved  glass-epoxy  panels  under 
load  control. 

Objective  of  Present  Work 

The  study  discussed  in  this  report  is  intended  to  extend  the  pres¬ 
ent  knowledge  of  the  static,  large  deflection,  behavior  of  orthotropic 
shallow  curved  panels.  To  accomplish  this,  an  exact  solution  is  presen¬ 
ted  for  a  simply  supported  cylindrical  panel  with  a  circular  cross 
section  and  subjected  to  a  line  loading  along  a  cylinder  generator.  The 
orthotropic  axes  are  parallel  to  the  generator  and  circumferential 
directions.  Derivation  of  the  equations  is  presented  in  Chapter  2.  Tne 
governing  equations  are  derived  using  laminated  plate  theory,  nonlinear 
strain-displacement  relations,  and  applying  variational  principles. 
Eccentricities,  or  offsets,  in  load  position  relative  to  midspan  are 
taken  into  account.  Due  to  the  geometry  of  the  problem  under  considera¬ 
tion,  only  a  one-dimensional  analysis  is  developed.  The  one  independent 
variable  is  related  to  the  circumferential  arc  length  along  the  panel . 
Because  of  the  one-dimensional  nature  of  the  problem,  ordinary  differen¬ 
tial  equations  govern  the  behavior  of  the  panel.  These  equations  are 
solved  in  closed  form.  The  aim  of  this  work  is  to  investigate  the 
effects  of  orthotropic  material  properties  on  the  behavior  of  a 
shallow  curved  panel.  In  past  work,  as  has  been  discussed,  the  behav¬ 
ior  of  the  structure  is  assumed  to  be  controlled  solely  by  the  geome¬ 
try.  For  an  orthotropic  curved  panel,  the  material  properties  may  also 
affect  the  response.  This  idea  is  investigated  here.  In  addition,  a 
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load-position  sensitivity  analysis  is  presented  along  with  predictions 
of  load  vs.  displacement,  load  vs.  thrust,  and  load  vs.  surface 
strains.  These  predictions  are  presented  in  Chapter  3.  In  order  to 
verify  the  analysis,  experimental  results  of  displacement-controlled 
tests  performed  on  graphite-epoxy  curved  panels  are  presented.  Chapter 
4  describes  the  experimental  set-up  while  Chapter  5  presents  experimen¬ 
tal  results.  Since  the  tests  are  run  under  displacement  control,  a 
greater  portion  of  the  equilibrium  path  is  stable  and  predictions  of  the 
buckling  load  under  load  control  may  be  made  with  greater  confidence. 
Also,  experimental  observation  of  the  post-buckling  mode,  unstable  under 
load  control,  may  be  made  using  displacement  control.  By  comparing 
analytic  and  experimental  results,  the  effectiveness  of  laminated  plate 
theory,  in  combination  with  an  analysis  accounting  for  large  elastic 
deformations,  in  predicting  the  structural  response  will  be  shown.  In 
addition,  the  influence  of  panel  geometry  and  panel  material  properties 
on  the  response  will  be  illustrated.  Finally,  Chapter  6  presents  con¬ 
clusions  and  recommendations  for  future  work. 
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n 

able  0.  Thus  in  the  analysis,  vv  =  w(e)  and  u  =  u(C-i).  Other  assumptions 
to  be  used  in  the  analysis  include  the  assumptions  of  a  linear  elastic 
material  behavior  and  linear  strain  variation  through  the  thickness. 

Also,  the  panel  is  assumed  to  be  thin  with  the  transverse  normal  stress 
negligible  in  comparison  with  the  inplane  stresses.  Transverse  shear 
deformations  are  neglected. 

Kinematics  of  Deformation 

The  kinematics  of  the  middle  surface  deformation  are  shown  in  fig. 
4.  The  applied  load  causes  each  point  of  the  panel  middle  surface  to 
displace  by  an  amount  t  and  the  corresponding  cross  section  to  rotate  by 
an  amount  Q.  The  unit  vectors  to  be  used  for  the  undeformed  geometry 
are  i  .  normal  to  the  middle  surface  and  directed  toward  the  center  of 

p» 

curvature,  and  the  tangential  unit  vector  i^.  For  the  deformed  config- 

A  A. 

uration  the  radial  and  tangential  unit  vectors  are  e^  and  e^,-  respect¬ 
ively.  The  position  vector  of  a  point  on  the  undefonned  middle  surface 

is  given  by, 

ft  =  -  Ri^.  (1) 

The  position  vector  for  the  same  point  after  deformation  is  denoted  by, 

r  =  ft  +  L  (2) 

The  displacement  vector  is 

X  =  w(S)ip  +  u(S)ij.,  (3) 

where  the  undeformed  arc-length  is 

S  =  Re.  (4) 

The  length  of  a  differential  element  of  the  undeformed  arc  length 
is  denoted  by  dS.  The  length  of  this  differential  element  of  arc  length 
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f"  the  defonnM  configurotfon  Is  represented  by  dS*.  By  definition  tne 
extensional  strain  of  the  middle  surface  is, 

0  dS*-dS 

^6  dS  '•  (5) 


To  derive  the  xjoverning  equations  in  terms  of  displacements,  it  is 
necessary  to  find  the  relation  betv/een  the  arc  lengths  and  displace¬ 
ments.  Ihis  relation  can  be  deterinined  by  realizing. 


From  eq.  2  the  change  of  the  deformed  position  vector  with  respect, to 
the  undeformed  arc  length  is. 


dr  df? 
d?  =  + 


dX 

d^’ 


(7) 


where. 


d^  c 
dS  ~  ’t* 


(8) 


and , 


To  obtain  eqs.  8  and  9  use  was  made  of  the  relations. 
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W’  " 


1  r 

4  1.  and  -~ 
R  t  d:j 


(10).  (11) 


Realizing  that  +  -^  represents  the  initiaT  curvature,  the  displace¬ 
ment  gradients,  and  r^,  are  defined  as. 


^n  = 


dw 


and  r^. 


_  du 
=  d^  ■ 


KqW 


(U).  (13) 


so  that  eq.  7  becomes, 

^  ^  (1  +  ^t^’f 

From  eqs.  6  and  14  the  deformed  arc  length  is  related  to  the  undeformed 
arc  length  by  the  relation. 


dS*  =  [{1  +  r,.)^  +  (15) 

The  extensional  strain,  e°,  is  now  calculated  to  be, 

O 

=  [(1  +  r^.)^  +  1  .  (16) 

Squaring  the  first  term  in  the  radical  and  applying  the  binomial  expan¬ 
sion  leads  to, 

"e  =  ^t  ^  +  o(r2). 


(17) 


By  assuming  tne  displacement  gradients  remain  small,  i.e., 

0  <  |r,.l  «  1  and  0  <  |rj  «  1,  (18).  (19) 

terms  of  third  order  and  higher  are  dropped  from  the  strain  expres¬ 
sion.  This  results  in  the  approximate  strain-displacement  relationship, 

0  du  1  ,dw  ,2 

^0  ~  dS  '  '^o”  I  ^dS  •  (20) 

The  configuration  of  the  deformed  middle  surface  is  specified  by 
the  position  vector  r(S*).  Let  <($*)  represent  the  curvature  of  the 
normal  section  of  this  surface  perpendicular  to  the  x  axis.  Then  the 
differential  geometry  of  this  normal  section  requires 


dr 

- 

dS*  " 

"t 

(21) 

A 

de 

p 

(22) 

A 

'‘‘'t 

•cer 

(23) 

These  expressions  are  analogous  to  eqs.  8.  10,  and  11  fortle  undeformed 
normal  section  of  the  middle  surface.  The  unit  tangent  vector,  eq.  21, 
is  related  to  the  displacement  gradients,  eqs.  12  and  13,  by  using  the 
chain  rule  to  relate  S*  to  S,  and  then  using  eqs.  5  and  14.  The  result 
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=  [I'nV  ^  • 

(24) 

Now  define  the  clockwise  angle  of  rotation  of  the  normal 
deformation  as  Q;  see  fig.  4.  Thus 

section  due  to 

=  sin (q) Ip  +  cos(Q)i^  » 

(25) 

and 

=  cos(Q)ip  -  sin(fi)i^  . 

(26) 

Comparing  eqs.  24  and  25,  the  trignometric  functions  of 

the  rotation 

angle  are 

sin(Q)  =  — ^ 

Uel 

(27) 

cos (Q)  =  - Q 

(2^) 

tdn(Q)  = 

(29) 

Consider  the  deformation  of  a  surface  parallel  to  the  middle  sur¬ 
face  specified  by  a  nonzero  constant  value  of  z.  Let  dS^  and  dS^* 
designate  the  undeformed  and  deformed  differential  arc-lengths,  respect 
ively,  of  the  normal  section  to  the  parallel  surface.  Consistent  with 
the  definition  of  the  middle  surface  strain,  eq.  5.  the  circumferential 

strain  of  the  parallel  surface  is  defined  as> 

• 
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^9  "" 


dS^  -  dS^ 


(30) 


dS, 


A  material  point  located  by  the  position  vector  on  the  undeforined 
parallel  surface  is  located  by  position  vector  on  the  deformed  par¬ 
allel  surface.  It  is  assumed  that  normals  to  the  undeforined  middle 
surface  remain  normal  and  undeforined  with  respect  to  the  middle  surface 
during  deformation  (Love-Ki rchhoff  hypothesis).  Thus  these  position 
vectors  are 

=  -  {R+z)i^,  and  =  r  -  ze^  .  (31), (32) 

Using  eqs.  8,  10,  21,  and  22,  the  differentials  of  these  vectors  with  z 

1 

held  constant  are 

dlT^  =  -  (i+ZKjj)dS  i^.,  and  dr^^  =  -  (l+ZK)dS*  e^.  (33), (.34) 

The  magnitudes  of  these  differential  vectors  are 


dS  =  (l+zic„)dS,  and  dS*  =  (l+ZK)dS* 


(35), (36) 


Substitute  these  into  eq.  30,  and  use  eq.  5,  to  get 


1+Z<  1 


(37) 


For  thin  shells, 
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0  <  IzK^I  <  IhK^I  <  <  I  .  (38) 

Assume  the  shell  is  thin,  and  expand  the  denominator  of  eq.  37  in  a 
power  series  in  z  using  the  binomial  series  to  get 

eg  =  +  0(z2)]  .  (39) 

Consequently,  for  thin  shells  and  small  strains  the  parallel  surface 
strain  is  approximated  by. 


eg  =  Eg  +  z(<-K^)  .  (40) 

This  result  shows  the  parallel  surface  strain  is  determined  by  the 
middle  surface  strain  and  its  change  in  Curvature. 

To  complete  the  kinematics,  the  change  in  curvature  of  the  middle 
surface  is  needed.  According  to  eq.  23  the  curvature  of  the  normal 
section  of  the  deformed  middle  surface  is 


K 


V'  dS*  * 


(41) 


Using  the  chain  rule  and  eq.  5  this  becomes 


=  e. 


1 


K 


(42) 
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From  eqs.  10,  11,  25,  and  26,  the  derivative  of  the  deformed  tangent 
unit  vector  with  respect  to  the  undeformed  arc  length  is. 


The  implication  of  the  approximation  In  eq.  46  is  important.  If  the 
rotations  are  assumed  to  be  small,  and  if  the  load-deformation  behavior* 
of  a  cur.ed  panel  from  no-load  to  complete  snap-through  is  to  be  invest¬ 
igated,  then  the  initial  curvature  of  the  panel  must  be  small.  This 
implies  the  work  here  is  limited  to  shallow  shell  theory.  The  curva¬ 
ture,  eq.  44,  can  now  be  expressed  as 


/ 

/  x- 
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"o  Hr  ’ 


(47) 


noting  that  the  small  strain  assumption  leads  to  the  approximation. 


1 


(48) 


Substituting  for  the  displacement  gradient,  the  curvature-displacement 
relationship  can  be  determined  from  eq.  47  as 


K 


(49) 


From  shallow  shell  theory  the  tangential  displacement  u  is  assumed 
to  be  small  relative  to  so  that  the  normal  displacement  gradient 
becomes, 

dw 

^n  =  H?*  (50) 


Thus  for  shallow  panels,  the  strain-displacement  and  curvature-displace¬ 
ment  relations  are. 


0  du 

^0  "  HF  ■ 


,dw.2 


(51) 


K 


.  d^w 
— 2 


(52) 


The  change  in  middle  surface  curvature,  K,  is  defined  to  be 

K  =  K  -  K^=  d^w/dS^ 


(53) 
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Total  Potential  Energy 

As  previously  stated,  the  governing  equations  will  be  derived  using 
variational  principles.  To  do  so  it  is  necessary  to  form  the  total 
potential  energy  of  the  system.  By  definition  the  total  potential 
energy  is, 

n  =  U  +  V  (54) 


where  U  is  the  strain  energy  stored  in  the  structure  and  V  is  the  poten¬ 
tial  of  the  loading.  The  strain  energy  density,  U(j,  is  defined  as. 


U  .  =  -:;r  Cf.  .  e.  .  , 
d  2  IJ  IJ 


Since  transverse  shear  strains  and  the  transverse  normal  stress  are 
assumed  to  be  negligible  in  comparison  to  the  inplane  stresses  and 
strains,  the  strain  energy  density  reduces  to, 


X  X  x0'xe  e 


The  previously  mentioned  one-dimensional  nature  of  the  geometry  suggests 
that, 

=  0  and  u^  =  0,  (57) 


where  is  the  displacement  in  the  x  direction.  Using  these  assump¬ 
tions  the  strain  energy  density  is  reduced  further  to. 
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U  =  -  a  E  .  (5«) 

2  0  0 

Since  the  varidbles  are  independent  of  x,  the  total  strain  energy,  using 
eq.  58,  is, 

•j  =  -r 

where  is  the  length  of  the  panel  in  the  x  direction. 

The  Constitutive  law  for  an  orthotropic  material  in  the  state  of 
plane  stress  described  here  is. 


(60) 


where  the  TJ  matrix  is  ttie  transformed  reduced  stiffness  matrix.  Equa¬ 


tion  60  is  written  in  structural  coordinates.  For  a  more  complete 
treatment  of  orthotropic  constitutive  relations  see  [9].  The  stress  in 
the  9  di rection  is. 


°Q  ^22^0* 

Substituting  this  relation  into  eq.  59  yields. 


(61) 


ORi’Gfr.'AL  P^.rz  f3 
OF  POOR  QUALilT 


/  / 

S  z 


Q22e0dSdz 


Strain  away  from  tiie  middle  surface,  from  eq.  40,  is  written  i 
the  middle  surface  strain  and  curvature  change  as, 


With  this  expression  the  constitutive  law  is, 

Integrating  this  expression  with  respect  to  z  yields, 


^22^0 


+  B22K 


the  force  resultant,  Ng,  being  defined  as. 


Ng  =  ogdz  . 


The  extepsional  stiffness,  A22,  and  bending-stretching  coupling 
ness,  B22»  are  defined  by. 


^22  “  ^22^k^^k  ■  ^k-1^ 


*^22  "  I  1^^^  ^^22^k  ^^k  ■  ^Ll^  ’ 


(52) 

terns  of 

(63) 

(64) 

(65) 

(66) 

stiff- 

(67) 

(68) 
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where  k  is  the  lamina  number  and  Zj^  and  (zj,  >  Z;.^i)  define  the 

location  of  the  kth  lamina  in  the  laminate.  Multiplying  the  constitu¬ 
tive  law  of  eq.  64  by  z  and  integrating  over  z  yield  the  expression  for 
the  moment  resultant.  Mg, 


where 


*^0  "  ®22^e  ^22*^’ 


Mg  =  OgZdZ. 


(69) 

(70) 


The  bending  stiffness,  defined  as, 

^22  ^  T  ^^22^k^^k  ■  ^k-1^  * 

Since  only  symmetric  lay-ups  are  to  be  examined,  the  stiffness  B22  is 
set  equal  to  zero  in  the  analysis. 

Substituting  the  strain  expression  of  eq.  63  into  eq.  62  and  inte 
grating  over  z  yields. 


Taking  the  first  variation  of  the  strain  energy  and  substituting  the 
constitutive  relations  of  eqs.  65  and  69  into  eq.  72  leads  to. 


6U  =  L  / 

S 


{Ng6e°  +  Hg6K}dS  . 


(73) 
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ORiG^fiA'.  r-r:; 
OF  POOR  QUA 


(82) 

(83) 

either  Mq  =  0  or  ^  Is  known.  (84) 

The  transition  conditions  at  S  =  » 


either  Ng  =  0  or  u  is  knovm, 


either  Ng  “  0  or  w  is  known. 


u  is  continuous,  (85) 

w  is  continuous,  (86) 

is  continuous,  (87) 

Mq  is  continuous,  (88) 

N  is  continuous,  (89) 

0 


dMg(^)  dMg(T-) 

— dS  (JS 


(90) 


Equation  90  indicates  that  the  shear  force  experiences  a  jump  at  S  =  S. 

An  additional  equation  can  be  obtained  by  integrating  both  sides  of 
eq.  51  over  the  arc  length  <  S  <  S^j,  Using  eq.  65,  recalling  822  = 

0,  eq.  51  assumes  the  form 


1  rdw->2 

7  • 


(91) 
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From  eq.  80,  N.  is  constant  with  S.  Therefore  integration  of 

.  9 

eq.  91 

results  In 

„  ^22  r^o  r  .  ^  1  fdW',2i -c 

(92) 
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• 

Where  u  vanishes  at  the  fixed  ends.  Using  eq.  80  in  eq.  81, 

results  in 

«  » A/\.„ 

(93) 

Equations  92  and  93  are  the  equations  governing  the  deformation  of 

the  panel.  For  the  panel  configuration  in  question  (see  fig. 

3)  the 

boundary  conditions  are. 

u  =  0 

(94) 

w  =  0 

(95) 

O 

n 

CD 

(96) 

at  S  =  ±  Sq.  Transition  conditions,  applied  at  S  =  S,  are. 

u(r)  =  u(S^) 

(97) 

[ 

wCT)  =  w(?") 

(98) 

.  . ■" 

(99) 

M0(S-)  = 

(100) 

NqCT)  =  Ngir") 

(101) 

dM  (T)  dM  (?") 
dS  dS  ^  • 

(102) 
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where  x  is  a  nondiinensional  parameter  defined  as. 


2  4 


(112) 


Equation  110  corresponds  to  eq,  93  while  eq.  Ill  corresponds  to  eq.  92. 

By  using  the  definitions  of  the  bending  and  extensional  stiffnesses 
for  an  isotropic  material,  the  parameter  \  is. 


(113) 


It  is  seen  that  the  material  dependence  is  eliminated  and  X  is  a  func¬ 
tion  of  geometry  only.  In  contrast,  with  an  orthotropic  material  X  is 
also  a  function  of  the  material  properties.  It  should  be  emphasized 
that  the  solution  to  eqs.  110  and  111  depend  only  on  the  parameter  X. 
since,  for  an  orthotropic  material,  X  incorporates  both  material  and 
geometric  properties  of  the  panel,  all  panel  characteristics  influence 
this  single  parameter. 

The  nondiinensional  boundary  conditions  and  transition  conditions, 
using  eqs.  105-109,  are. 


u  =  0  at  r  =  ± 

1 

(114) 

(t)  =  0  at  r  ± 

1 

(115) 

=  0  at  r  = 
dr 

±  1 

(116) 

u  o 
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u(r)  =  uo^) 


<!>(?■)  =  (i>(r^) 
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nCT)  =  r,(r) 


d-tO  =  .^■-».(pL  + 


(118) 


(119) 
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Chapter  3 

SOLUTION  TO  GOVERNING  EQUATION  AND  NUMERICAL  RESULTS 
The  governing  differential  equation  developed  in  the  previous 
section,  eq.  110,  appears  to  be  a  linear  ordinary  differential  equation 
with  n  being  an  unknown  constant.  This  is  not  the  case  since  q  is 
actually  a  nonlinear  function  of  the  displacement.  However,  linear 
analysis  may  be  used  to  find  the  functional  form  of  the  radial 
displacement,  <)(r).  The  nonlinearity  of  the  problem  manifests  itself 
through  the  second  governing  equation,  eq.  111.  Since  the  functional 
form  of  (t)(r)  is  established  by  eq.  110,  eq.  Ill  results  in  a  transcen¬ 
dental  equation  for  q.  Since  only  compressive  middle  surface  strains 
(i.e.,  Nq  <  0)  are  of  interest,  the  thrust  parameter  q  is  defined  such 
that. 


q^  =  -  q\^  (123) 

Substituting  this  expression  into  the  governing  equations  yields, 

A  *  *  1)  .  0  (124) 

dr  dr 

Equation  124  can  be  solved  in  closed  form.  However,  the  form  of  the 
solution  depends  on  the  value  of  q.  The  solutions  for  various  values  of 
q  are  discussed  below. 


32 


33 


Midspan  loading:  q  /  O  and  g  /  n  =  1,  3.  5,  ... 

With  the  load  acting  radially  at  midspan,  T  =  0,and  with  q  /  0  and 

q  where  n  =  1,  3,  5,  ....  the  general  solution  to  eq.  124  for  the 

radial  displacement  is 


({)(r)  =  +  A^r  -  +  a^  cos(qr)  +  a^  sin(qr)  (126) 

In  each  of  the  intervals  -1  <  r  <  0  and  0  <  r  <  1.  Applying  the  bound¬ 
ary  conditions,  eqs,  114-116,  on  each  interval,  and  the  transition 
conditions,  eqs.  117-122,  at  F  =  0  yields. 


♦  (r)  .  b  sln(qr)  ♦  ^2  (1  -  -—If})  +  I  (1  -  r^) 

4 

-  [1  -  |r|  -  i-  tan(q)cos(qr)  +  -i  sin(q|r|)}  . 


(127) 


This  solution  is  valid  over  both  intervals.  Iri  this  expression  b  repre¬ 
sents  the  magnitude  of  the  asymmetric  displacement  response.  The  fol¬ 
lowing  conditions  on  b  are  determined  from  the  boundary  condition 

i 

requiring  the  moment  to  vanish  at  the  edges: 

if  q  nit  then  b  =  0 
if  q  =  nil  then  b  /  0, 

where  n  =1,  2,  3 .  To  complete  the  solution  process  the  func¬ 

tional  form  of  ({)(r),  eq.  127,  is  substituted  into  eq.  125.  This  relates 
the  thrust  parameter  q  to  the  nondimensional  applied  load,  p.  After 
some  algebra  this  relationship  is  written  as. 
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+  C  =  0  » 

where  A,  B,  and  C  are  functions  of  q  and  K,  and 

e  = 

Appendix  A  contains  the  functional  form  of  A,  B,  and  C  for  all  cases 
considered  here.  For  a  given  panel ,  \  is  specified.  In  eq.  128,  by 
specifying  the  thrust  parameter  q,  a  quadratic  equation  can  be  solved  to 
obtain  the  corresponding  load,  p.  This  is  in  opposition  to  solving  eq. 
128  as  a  transcendental  equation  for  the  thrust  parameter  if  loads  are 
first  specified.  In  this  solution  procedure  it  is  important  to  recog¬ 
nize  that  eq.  128  may  have  two  real  roots,  one  real  root,  or  no  real 
roots,  depending  on  the  values  of  \  and  q.  Each  real  root  corresponds 
to  a  separate  equilibrium  configuration.  Having  found  the  relationship 
between  load  and  thrust,  the  load-displacement  relationship  follows  from 
eq.  127. 

For  q  =  nit,  b  0,  there  is  an  equation  similar  to  eq.  128  which 
relates  the  applied  load  to  the  magnitude  of  the  asymmetric  deformation, 
b.  It  also  originates  from  eq.  125.  This  relationship  is  given  by  the 
expression, 

l)b^  :=  +  B^  +  C.  (130) 

Here  A,  B,  and  D  are  constants  depending  only  on  the  value  of  n.  Con¬ 
stant  C  depends  on  n  and  X.  Again  these  constants  are  expressed  in 
Appendix  A.  Real  solutions  for  the  asymmetric  amplitude  b  require  the 
right-hand-side  of  eq.  130  to  be  nonnegative.  Consequently,  not  all 


(128) 

(129) 
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integer  values  of  n  are  permissible  for  a  specific  value  of  Notice 
that  two  asymmetric  configurations  are  determined  from  eq.  130  at  given 
load  if  the  right-hand-side  is  positive. 

Hidspan  Loading:  q  =  0 

When  the  thrust  is  equal  to  zero  eq.  127  is  not  the  correct  func¬ 
tional  form  of  4)(r).  Setting  q  =  0  in  the  governing  equation,  eq.  110, 
a  new  functional  form  of  (t)(r)  is  obtained.  Applying  the  boundary  condi-i 
tions  (eqs.  114-116)  and  the  transition  conditions  (eqs.  117-122)  the 
functional  form  of  the  radial  displacement  is, 

<i>(r)  =  I  "  "I"  ■*■  ■§■  }  •  (131) 

By  substituting  eq.  131  into  eq.  125  the  relationship, 

2  25 

=  (132) 

is  obtained  for  determining  the  loads  at  which  the  thrust  parameter,  and 
thus  the  middle  surface  thrust,  N^,  vanishes.  This  occurs  at  either 

D 

25 

C  =  0  or  C  =  -^  .  The  C  =  0  solution  states,  as  expected,  that  there 
is  no  thrust  in  the  panel  if  there  is  no  load.  However,  the  nonzero 
solution  to  eq.  132  indicates  that  there  is  a  loaded  and  deformed  state 
of  the  panel  for  which  there  is  no  thrust  in  the  panel.  In  this  special 
situation  the  load  is  resisted  by  bending  only. 
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lidspan  Loading:  =  gl,  „  =  1.  3,  5, 

The  functional  form  for  (i,(r)  in  eq.  127  is  incorrect  when  q  =  "i 
n  =  1,  3,  5,  ,  In  this  case  the  linear  inhomogeneous  mathematical 

problem  for  $(r)  ,  given  by  eq.  124,  the  boundary  conditions,  and  the 
transition  conditions,  does  not  have  a  unique  solution.  The  mathemati¬ 
cal  reason  for  this  is  found  in  Fredholm's  alternative  theorem  [10,  pg. 
547].  The  associated  adjoint  homogeneous  problem  for  (i)(r)  has  nontriv¬ 
ial  solutions  (eigenfunctions)  associated  with  these  values  of  q  (eigen¬ 
values).  Consequently  the  tonditions  for  a  solution  to  (|>(r)  require, 

+1  2 

c  <l>^(0)  =  /  ^q  (!.^(r)dr  ,  (133) 


where  4>|^(r)  =  cos(^H-),  m  =  1,  3,  5,  ,..  are  the  eigenfunctions  of 

the  associated  homogeneous  problem.  (Equation  133  is  equivalent  to  tne 
"inhomogeneity"  being  orthogonal  to  the  eigenfunctions.)  Equation  133 
determines  the  load  for  the  given  values  of  the  thrust.  The  displace- 
merit  functions  which  result  from  this  are. 


,  V  1  2  *  2sin(^) 

»'n(^)  =  2  (1  -  r  )  +  cos(qr),+  — — (|rl  -  i) 


riTc 


sin(q|r|)}  , 


(134) 


n  =  1,  3,  5,  ...,  n^^jj  .  The  values  of  A*  are  found  by  substituting  eq. 
134  into  eq.  125.  This  yields 
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A(A*)^  +  BA*  +  C  =  0  . 


(135) 


where  A,  B,  and  C  are  constants  depending  on  the  values  of  n  and  x. 
Appendix  A  presents  the  coefficients  A,  B,  C.  The  solution  for  ^  is  not 
unique  In  this  case.  The  finite  number  of  these  solutions, 
follows  from  the  requirement  of  real  solutions  to  eq.  135. 

Eccentric  Loading:  g  7^  0  and  q  n  1,  2,  3.  ... 

For  loads  acting  radially  at  the  location  r  =  T,  0,  the  general 

solution  for  the  radial  displacement,  with  q  /  0,  q  /e  is 


<t.(r)  =  ^  (1 


COS 


cos(q 


)  +  -i  (1  -  r^)  -  {1  +  r  sgn(r) 


2q 


-  r  r  -  |r|  -  ~  tan(q)  cos(qr}  cos(qr) 


+  ^sgn(r)  sin(qr)  cos(qr)  -.^^sin(qr)  sin(qr) 


(136) 


where 


+  -  cos(qr)  Sin(q|r|)  , 


-1,  r  <  0 
sgn(r)  =  0,  r  =  0 

1,  r  >  0  . 


(137) 


Equation  136  satisfies  the  boundary  conditions  and  the  transition  condi¬ 
tions.  Substitution  of  eq.  136  into  eq.  125  yields  the  load-thrust 
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relationship, 

+  C  =  0 


(138) 


where  A,  B,  and  C  are  functions  of  the  thrust,  k,  and  r.  (See  Appendix 
A).  The  solution  procedure  for  the  offset  loading  follows  that  of  the 
inidspan  loading  case. 


Eccentric  Loading:  g  =  0 

For  the  thrust  parameter  equal  to  zero,  the  governing  equation  is 
again  altered.  The  resulting  displacement  solution  is. 


4>(r)  =  (r  +  sgh(r))  -  ^  (l  +  r  sgn(r)) 


(139) 


The  loads  for  which  the  thrust  is  zero  are  given  by  the  relationship, 

5  +  l8{r  )^+  (T  f 

-  4  t - 4...;^]}  =  0  ,  (140) 

^  3-7(r  )^+  5(r  )  -{r  )'" 


Which  is  obtained  by  substituting  the  functional  form  of  the  radial 
displacement  into  eq*  125* 


Eccentric  Loading,  q  =  n  =  1,  2,  3> 

With  q  -^9  n  =  1,  2,  3,  ...  the  adjoint  problem  outlined  for  the 
midspan  loading  must  again  be  examined.  The  orthogonal ity  condition 
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A(An)^  +  BA^  +  C  =  0  ,  (I44) 

where  A,  B,  and  C  are  constants  which  depend  on  the  values  of  tne 
thrust.  X,  and  r  (See  Appendix  A  for  the  functional  forms  of  A,  B.  and 
C). 


Typical  Numerical  Results  -  Midspan  Loading 

ShaUow  cylindrically  curved  panels  exhibit  a  number  of  different 
deformation  and  buckling  behaviors.  The  particular  behavior  depends  on 
the  parameter  X.  Some  of  the  effects  of  orthotropic  material  properties 
on  the  value  of  X,  and  thus  the  behavior  of  curved  panels,  is  explored 
in  Appendix  B.  Generally,  though,  curved  panels  which  have  no  eccen¬ 
tricity  in  the  loading  (r  =  0)  exhibit  three  distinct  behaviors.  In 
this  section  each  type  of  behavior  will  be  discussed  in  order  to 
familiarize  the  reader  with  typical  analytic  results.  These  discussions 
will  focus  on  the  behavior  under  displacement  control  and  load  con¬ 
trol.  However,  the  condition  represented  by  the  experimental  setup  to 

be  described  is  displacement  control. 

If  X  is  sufficiently  small  the  load-displacement  behavior  is  mono¬ 
tonic  and  no  snap-buckling  occurs.  In  fig.  5  typical  load-thrust,  p  vs. 
q.  and  load-mi dspan  displacement,  p  vs.  ?,  diagrams  are  shown  for  a 
panel  i^hich  exhibits  monotonic  behavior.  Here  $  denotes  the 
nondiinensional  radial  displacement  at  the  location  of  the  applied 
load.  As  the  load  and  displacement  increase  from  zero  through  points  1 
and  2,  the  stiffness  (dp/d?)  of  the  panel  decreases  while  the  thrust 
increases.  At  point  3  the  thrust  is  a  maximum  and  the  stiffness  is  a 
minimum.  A  further  increase  i-  load  or  displacement  through  points  4 
and  5,  depending  on  which  variable  is  controlled,  causes  the  thrust  to 
decrease  and  the  stiffness  to  increase.  Finally  at  point  6,  a  point 
close  to  the  fully  inverted  configuration,  the  thrust  is  approaching 
zero.  In  fig.  6  the  equilibrium  configurations  corresponding  to  the 
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points  denoted  on  the  load-thrust  and  load-displacenient  are  shown. 

Figure  7  shows  the  behavior  of  the  upper  and  lower  surface  strains, 
i.e.,  z  =  +  Y  ^hd  ^  -  -j*  St  the  two  guarter  points  (r  =  +  0.5)  and 

midspan  (r  =  0,0)  ,  as  a  function  of  the  applied  load.  Obviously  due  to 
symmetry  of  the  response,  the  strains  at  the  quarter  points  are 
Identical,  This  panel  exhibits  no  instability,  a  characteristic  which 
is  typical  of  panels  with  \  <  1.9762. 

A  second  type  of  behavior  is  known  as  limit  point  behavior.  This 
is  depicted  in  fig.  8.  As  the  load,  or  displacement,  is  increased  from 
zero  through  point  1,  the  thrust  increases  and  the  stiffness  decrea-* 
ses.  At  point  2  the  load-thrust  curve  has  a  horizontal  tangent  and  the 
stiffness  of  the  panel  is  zero.  At  point  2a  relative  maximum  load,  or 
limit  point,  is  reached  on  the  equilibrium  path.  Under  load  control 
this  is  a  critical  equilibrium  configuration  and  increasing  the  load 
results  in  a  sudden  dynamic  snap-through  to  a  equilibrium  configuration 
between  points  5  and  6.  Under  displacement  control  this  point  is  sta¬ 
ble.  Increasing  the  displacement  past  point  2  will  cause  the  load  to 
decrease  and  the  thrust  to  continue  to  increase.  Point  3  denotes  the 
maximum  thrust  for  the  panel.  Increasing  the  displacement  from  point  3 
to  point  4  causes  both  the  thrust  and  the  load  to  decrease.  Notice  that 
point  4  is  also  a  load  limit  point.  From  point  4  through  point  5  the 
load  increases  while  the  thrust  decreases  to  zero  with  increasing  dis¬ 
placement,  Figure  9  shows  the  equilibrium  configurations,  under  dis¬ 
placement  control,  for  various  values  of  load  and  thrust.  Figure  10 
illustrates  the  strain  behavior  at  the  two  quarter  points  and  at  mid- 
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span.  Under  displacement  control  this  panel  exhibits  no  instability. 
Under  load  control  there  is  the  snap-through  instability  at  point  2,  the 
limit  pointy  which  is  considered  the  buckling  load.  Limit  point  behav¬ 
ior  is  found  for  1.9762  <  X  <  2.825. 

The  final  type  of  behavior  possible  for  the  panel  with  no  load 
eccentricity  is  known  as  bifurcation  behavior.  Bifurcation  occurs  on 
the  initial  loading  path  when  the  maximum  thrust  equals  r..  Figure  11 
shows  typical  load-thrust  and  load-displacement  diagram,  for  bifurcation 
behavior.  Load,  displacement,  and  thrust  all  increase  f  om  no-load 
through  point  1  to  point  2.  At  point  2  a  secondary  equr  orium  path, 
the  straight  line  through  points  2,  3,  and  4,  intersects  the  primary 
path.  This  secondary  path  represents  the  asymmetric,  response  mentioned 
in  Chapter  2  and  included  in  eq.  127.  (Actually  there  are  two  equilib¬ 
rium  paths  associated  with  asymmetric  configurations  intersecting  at 
point  2  since  two  values  of  the  amplitude  eq.  130,  are  determined  for. 
each  value  of  the  load.)  Point  2  is  called  a  bifurcation  point.  In 
load  control  both  the  primary  equilibrium  path,  beyond  point  2,  and  the 

I 

secondary  equilibrium  path  configurations  are  unstable.  Thus  the  panel 
will  dynamically  snap-through  at  the  bifurcation  point  to  a  configura¬ 
tion  on  the  primary  path  beyond  point  6.  Point  2  is  considered  the 
buckling  load.  In  displacement  control  the  secondary  path  is  stable 
while  the  primar.'  path  becomes  unstable.  Increasing  the  displacement 
causes  the  load  to  decrease,  along  the  straight  line  between  points  2, 

3,  and  4,  and  the  thrust  to  remain  constant.  Note  the  load  actually 
drops  to  zero  at  point  3.  The  spatial  configuration  of  the  panel 
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corresponding  to  points  along  the  secondary  path  is  asymmetric  while  the 
primary  path  corresponds  to  symmetric  deformation.  Point  4  represents 
the  second  bifurcation  point  where  the  secondary  path  again  intersects 
the  primary  path.  Figure  12  illustrates  the  spatial  deformation  charac¬ 
teristics  for  a  panel  which  exhibits  bifurcation  behavior  and  is  under 
displacement  control.  Note  the  asymmetry  in  the  deformed  shape  as  the 
response  moves  along  the  bifurcated  path.  The  surface  strains  shown  in 
fig.  13  illustrate  only  those  strains  predicted  for  the  stable  equi¬ 
librium  path  under  displacement  control. 

Bifurcated  solutions  exist  at  all  values  of  q  =  nit,  n  =  1,  2,  3, 

. ..  but  only  the  first  bifurcated  path  (n  =  1)  represents  the 

physically  important  result  for  the  static  response.  It  should  be  noted 
here  that  bifurcation  behavior  is  a  mathematical  result  of  an  analyses 
for  a  perfectly  symmetrical  structure.  Since,  in  practice,  real  struc¬ 
tures  are  not  perfect  this  idealization  will  never  be  achieved.  How¬ 
ever,  with  relatively  small  initial  imperfections  it  cbuld  be  closely 
approximated.  Bifurcation  behavior  occurs  for  \  >  2.825. 


Typical  Numerical  Results  -  Eccentric  Loading 

When  offsets,  or  eccentricities,  in  the  location  of  the  applied 
loading  are  present,  only  two  types  of  behavior  exist.;  These  are  mono¬ 
tonic  behavior  and  limit  point  behavior.  Panels  which  exhibit  monotonic 
or  limit  point  behaviors  when  the  load  is  at  midspan  show  little  change 
in  their  behavior  for  small  eccentricities  in  the  loading.  Figures  14- 
16  show  the  monotonic  response  of  a  panel  with  geometric  and  material 
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properties  identical  to  the  panel  of  figs.  5-7,  with  the  exception  of  a 
10%  offset  in  the  load  being  present.  Here  the  load  is  offset  to  the 
right.  The  only  significant  effect  of  the  offset  is  a  slight  asymmetry 
in  the  deformation.  This  is  illustrated  by  the  equilibrium  configura¬ 
tions  shown  In  fig.  15.  Figure  16  shows  the  surface  strain  behavior  for 
this  loading  case. 

Figures  17-19  show  the  effects  of  load  offset  on  a  panel  which 
exhibits  limit  point  behavior  when  the  load  is  exactly  at  midspan. 

Again  the  only  significant  effect  is  a  slight  asymmetry  in  the  deforma¬ 
tion.  A  comparison  of  figs.  8  and  17  shows  that  the  offset  leaves  the 
limit  load  essentially  unchanged, 

A  panel  which  exhibits  instability  at  bifurcation  points  when  the 
load  is  exactly  at  midspan  is  affected  to  a  much  greater  degree  by 
eccentricities  in  the  load  position.  This  is  easily  seen  by  comparing 
the  load-thrust  and  load-displacement  diagrams  of  fig.  11,  the  case  of 
midspan  loading,  with  those  of  fig.  20.  Following  the  equilibrium  path 
shown  in  fig.  20,  an  increase  in  the  displacement  causes  the  thrust  and 
load  to  increase  up  to  point  2.  Here  point  2  represents  a  load  limit 
point  (horizontal  tangent).  For  the  midspan  loading  case,  fig.  11, 
point  2  was  not  a  limit  point.  It  was  a  bifurcation  point.  For  the 
offset  loading,  any  increase  in  the  applied  load  beyond  point  2  would 
cause  loss  of  stability.  On  the  other  hand,  controlling  the  displace¬ 
ment  from  point  2  to  point  4  causes  the  load  to  decrease  and  the  thrust 
to  continue  increasing.  At  point  4,  q  >  it.  Thus  for  load  offsets  q 
=  It  is  no  longer  the  limiting  value  of  thrust  for  the  stable  equilibrium 


Note:  numbers  correspond  to  points  on  fi 
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path.  RecalJ  from  fig.  11,  q  =  -rt  was  the  limiting  value  for  the  midspan 
case.  This  allows  the  eccentrically  loaded  panel  to  exhibit  higher 
strain  levels  than  are  found  with  midspan  loading.  Also,  at  point  4  a 
displacement  limit  point  (vertical  tangent)  is  found.  Therefore  an 
increase  in  the  displacement  causes  the  arch  to  snap-through  to  point 
5.  The  load  level  p  suddenly  jumps,  as  does  the  thrust  level.  From 
point  5  an  increase  in  the  displacement  causes  the  load  to  increase 
while  the  thrust  decreases  to  zero.  A  comparison  of  figs.  20  and  11 
reveals  that  the  bifurcation  load  of  fig.  11  is  significantly  higher 
than  the  limit  load  of  the  offset  response  in  fig.  20.  Figure  21  snows 
the  equilibrium  configurations  for  various  points  along  the  equilibrium 
path.  Of  special  interest  are  the  shapes  associated  with  points  4  and 
5.  The  shape  corresponding  to  point  4  represents  the  equilibrium  con¬ 
figuration  just  before  the  displacement  limit  point  while  the  shape 
corresponding  to  point  5  is  the  configuration  of  the  panel  just  after 

I 

the  displacement  limit  point.  Figure  22  illustrates  the  behavior  of  the 
surface  strains  vs.  the  applied- load,  it  should  be  noted  here  that  the 
sharp  jumps  in  the  strain  curves  are  due  to  the  displacement  limit 

!  r 

point.  This  portion  of  the  curve  is  not  stable,  the  surface  strains 
experience  a  dynamic  change  due  to  the  change  in  Iqsd  and  thrust. 

From  the  analysis  of  the  response,  it  appears  that  the  behavior  of 
a  cylindrical  panel  is  quite  sensitive  to  load  eccentricities  for  cer¬ 
tain  ranges  of  X.  To  fully  understand  the  effect  of  load  offsets,  the 
relation  between  the  magnitude  of  the  buckling  load  and  the  offset  angle 
is  examined.  Since  the  problem  is  physically  symmetric  with  respect  to 


65 


tne  midspan  x'erticai,  only  positive  load  offsets  are  shown.  Figure  23a 
shows  the  effect  of  load  eccentricity  on  the  buckling  load  for  three 
different  panels,  all  of  v/hich  exhibit  limit  point  behavior  when  the 
load  is  exactly  at  midspan.  Here  the  limit  load  is  considered  the 
buckling  load.  For  panel  1,  which  has  X.  =  1.99,  increasing  the  offset 
angle  causes  only  slight  increases  in  the  buckling  load.  However,  with 
this  panel,  an  offset  causes  the  behavior  to  switch  from  limit  point  to 
monotonic  behavior.  Thus  there  no  longer  is  an  instability.  Panel  2, 
with  \  =  2.05,  resembles  panel  1  with  the  exception  that  limit  point 
behavior  applies  for  the  entire  range  of  load  offsets.  Panel  3, 
where  \  =  2,5,  is  different  than  the  first  two  in  that  the  buckling  load 
initially  decreases  with  increasing  load  offset.  Still  the  panel  is  not 
overly  sensitive  to  eccentricities  in  the  loading.  Small  perturba¬ 
tions  (f  «  0,1)  in  the  load  position  only  cause  negligible  changes  in 
the  buckling  load.  Both  panels  2  and  3  exhibit  buckling  loads  approach¬ 
ing  infinity  as  the  applied  load  is  moved  very  close  to  the  supports. 

A  panel  which  exhibits  bifurcation  instability  when  the  load  is 
exactly  at  midspan  tends  to  be  quite  sensitive  to  load  eccentricities. 
Figure  23b  shows  the  buckling  load  as  a  function  of  offset  angle  for 
such  a  panel.  The  bifurcation  load  corresponds  to  the  buckling  load 
for  T  =  0,  and  tne  limit  point  load  corresponds  to  the  buckling  load  for 
all  other  load  offsets.  The  buckling  load  shows  significant  reductions 
even  for  small  load  offsets.  This  type  of  behavior  leads  to  the  panel 
being  considered  imperfection-sensitive. 


Cnapter  4 

EXPEKIMENTAL  SET-UP  AND  PROCEDURE 

Little  experimental  work  has  been  performed  In  the  area  of  the 
large,  out-of-plane  deformation  behavior  of  shallow  composite  cylin¬ 
ders.  The  work  that  has  been  presented  tends  to  be  somewhat  limited  in 
scope  in  that  the  range  of  geometry  is  small  or  the  effects  of  load 
eccentric^ies  have  not  been  investigated.  To  check  the  validity  of  the 
analysis  just  presented,  a  series  of  tests  on  composite  specimens  oi 
varying  geometries  and  load  offsets  were  performed.  These  tests  were  to 
determine  if  all  the  behavioral  characteristics  predicted  by  the  analy¬ 
sis  are  observable.  Also,  the  effectiveness  of  the  solution  procedure, 
as  well  as  classical  lamination  theory,  in  determining  the  large  defor¬ 
mation  response  of  composite  structures  was  to  be  verified. 

In  order  to  make  a  comparison  of  the  analytic  and  experimental 
results,  the  conditions  of  the  experiment  needed  to  closely  approximate 
those  used  for  the  analysis.  This  required  the  test  specimen  to  be 
supported  such  that  the  ends  remained  free  of  rotational  constraints 
while  the  span  between  supports  Was  kept  constant.  Also  the  load  needed 
to  act  radially  at  the  panel  middle  surface.  To  meet  these  require¬ 
ments,  as  well  as  others,  a  test  fixture  and  test  procedure  were 
designed  to  facilitate  a  fair  comparison  of  theory  and  experiment. 

Curved  Panel  Specimens 

The  test  specimens  were  fabricated  by  NASA-Langley  using  AS4/3502 
graphite-epoxy  pre-preg.  Table  1  shows  the  specified  geometry  and 
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Table  1 


Intended  Test  Specimen  Geometry 


Specimen 

Number 

Lay  up 

Thickness 

(in-) 

Load 

Offset 

(%) 

1 

[(90/0)4]s 

0,08 

0 

2 

C(90/0)43s 

0.03 

15 

3 

C(90/0)4]s 

0.08 

30 

4 

[(90/0) 435 

0.08 

50 

5 

[(90/0)3], 

0.06 

0 

6 

[(90/0)3], 

0.06 

15 

7 

[(90/0)3], 

0.06 

30 

8 

[(90/0)3], 

0.06 

50 

9 

[90/45/0/-45], 

0.04 

15 

10 

[90/45/0/-45], 

0.04 

30 

11 

[90/45/0/-45], 

0.04 

50 

12 

[(90/0)4], 

0.08 

0 

13 

[(90/0)4], 

0.08 

15 

U 

[(90/0)4], 

0.08 

30 

15 

[(90/0)]^ 

0.02 

0 

16 

[(90/0)], 

0.02 

[5 

17 

[(90/0)], 

0.02 

30 

18 

[(90/0)], 

0.02 

50 

Note:  For  all 

cases  the  intended  radius 

v/as  60  in.  with 

an  arc  length 

of  12  in 

with  the  exception  of 

specimens  12-14, 

v/here  the  arc 

length  was  6  in. 
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stacking  sequence  for  each  specimen  tested.  For  each  geometry'  which  was 
to  be  tested,  a  large  18  x  24  in.  curved  panel  v/as  fabricated  using 
existing  curved  forms.  After  curing,  each  large  panel  was  ultrasonic- 
ally  C  -scanned  to  determine  if  any  irregularities  were  present.  For 
each  large  panel  the  results  of  the  C-scan  showed  uniformity  in  the 
transmittance  of  sound  waves,  suggesting  a  uniformity  of  the  material. 
Upon  determining  the  quality  of  the  panels,  the  specimens  were  then 
machined  from  the  larger  panels.  Each  specimen  had  a  set  of  three  holes 
drilled  at  the  point  of  load  application.  This  will  be  discussed 
shortly.  By  keeping  a  tight  tolerance  on  the  arc-wise  location  of  the 
holes,  it  was  possible  to  control  the  amount  of  load  eccentricity. 

After  machining,  slotted  dowels  0.375  in.  in  diameter,  to  be  further 
discussed  in  connection  with  the  test  fixture,  were  attached  to  each  end 
of  the  specimens.  Figure  24  shows  6  in,  arc  length  and  12  in.  arc 
length  panels  both  with  and  without  the  dowels  attached.  Two  dov/els  are 
shown  in  fig.  24,  as  are  two  larger  dowels  which  will  be  discussed. 
Strain  gaging  has  been  attached  to  the  panels  in  fig.  24  and  a  6  in, 
rule  is  visible  in  the  photograph. 

Test  Fixture 

Figure  25  snows  a  schematic  of  the  test  fixture  wnile  fig.  26  shows 
a  photograph  of  the  fixture.  The  major  components  of  the  fixture 
are:  a)  end  dowels,  b)  support  blocks,  c)  a  base  plate,  d)  loading 

head  and  pivot,  and  e)  tongue  and  connector.  Each  component  and  its 
function  will  now  be  described. 


FI6.  25  Schematic  OF  TEST  fixture 
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a)  End  Dowels 

The  end  dowels  were  0.375  in.  diameter  steel  rods,  slotted  to 
accommodate  the  ends  of  the  curved  panels.  The  dowels  were  one  part  of 
a  system  designed  to  accurately  represent  a  simple  support  condition. 
The  ends  of  the  panel  were  potted  solidly  in  the  dowel  slots  using 
epoxy.  The  epoxy  prevented  the  ends  of  the  panel  from  brooming  under 
the  inplane  compressive  thrust  loads.  This  potting  is  shown  in  fig. 

27.  The  rods  were  slotted  along  a  rod  radius  and  the  panels  were  Cen¬ 
tered,  using  shims  if  necessary,  within  the  slots.  With  care,  the  end 
of  the  middle  surface  of  the  panel  could  be  made  to  coincide  with  the 
centerline  of  the  dowel.  With  this  arrangement,  the  boundary  conditions 
at  the  ends  of  the  middle  surface  (S  =  ±Sq)  could  be  accurately 
enforced. 

b)  Support  Blocks 

The  support  blocks  were  an  important  component  in  the  system  and 
two  designs  were  used  before  the  test  fixture  was  satisfactory.  The 
support  blocks  provided  tiie  simple  support  condition  by  preventing 
translational  displacements  but  allowing  rotational  displacement.  The 
original  design  called  for  support  blocks  into  which  semicircular 
grooves  were  machined.  The  grooves  were  designed  to  accommodate  the 
circular  dowels.  A  lubricant  between  the  dowel  and  the  groove  was  used 
to  minimize  frictional  effects  related  to  rotation  of  the  ends  of  the 
specimen.  Unfortunately,  due  to  the  inplane  thrust  forces  induced  in 
the  curved  panel  as  it  deformed,  significant  frictional  forces  developed 


in  the  end  blocks.  These  frictional  forces  resulted  In  end  moments  on 
the  panel,  a  violation  of  the  moment-free  boundary  condition  assumed  in 
the  analysis.  A  second  design  of  the  end  block  used  six  roller  bearings 
in  each  block.  These  roller  bearings  were  situated  so  that  the  dowels 
were  forced  against  them.  This  is  shown  in  fig.  28.  There  were 
actually  three  bearings  in  a  line  above  the  end  dowel  and  three  in  a 
line  below  the  end  dowel.  Figure  29  shows  a  photograph  of  the  roller 
bearing  end  blocks. 

c)  Base  Plate 

Figure  30  shows  a  schematic  of  the  base  plate.  The  base  plate  is 
also  visible  in  figs.  25  and  26.  The  function  of  the  base  plate  was  to 
provide  stiffness  to  the  entire  fixture  and  to  provide  a  solid  founda¬ 
tion  for  locating  the  support  blocks.  The  base  plate  was  19  x  4  x  1  in., 
and  was  machined  from  a  larger  piece  of  steel  stock.  Care  was  taken  to 
insure  that  the  top  and  bottom  surfaces  of  the  plate  were  flat  and 
parallel  with  each  other.  One  support  block  was  mounted  at  the  location 
of  the  six  isolated  holes  (top  of  fig.  30)  while  the  other  was  mounted 
in  either  set  of  slots.  Hardened  steel  pins,  0.250  in.  in  diameter  and 
protruding  from  the  base  plate  into  the  support  block,  were  also  used  to 
secure  the  block  mounted  at  the  holes.  More  will  be  said  of  this 
shortly.  The  two  sets  of  slots  were  to  accommodate  the  panels  df  two 
different  lengths.  The  slots  were  machined  into  the  plate  so  that  the 
position  of  the  support  blocks  relative  to  one  another  could  be 
adjusted.  This  was  necessary  since  the  panels  all  had  slightly  differ- 
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FIG.  29  Photograph  of  support  block 


ent  lengths  when  machined.  The  block  at  the  slotted  end  could  be 
adjusted  along  the  length  of  the  base  plate  so  that  the  panel  was  tight 
in  the  fixture  but  without  the  radius  of  curvature  being  too  different 
than  the  natural  value,  i.e.,  the  value  specified  in  Table  1.  This  was 
very  important  since  the  response  of  curved  panels  can  be  quite 
sensitive  to  the  initial  radius  of  curvature. 

There  was  concern  that  the  large  inplane  forces  generated  as  the 
panels  defonned  would  cause  the  blocks  to  move  relative  to  the  base 
plate.  This  is  the  reason  hardened  steel  pins  were  used  with  the  one 
support  block.  For  the  support  block  attached  at  the  slots,  another 
design  was  initially  incorporated.  A  third  steel  block  was  bolted  into 
position  at  the  six  holes  near  the  slots.  This  block  was  designed  to  be 
a  back-up  support  for  the  block  at  the  slots.  The  load  was  transferred 

to  the  back-up  by  way  of  several  horizontally  oriented bolts  jamnied 

I 

tjetween  the  back-up  block  and  the  support  block.  As  it  turned  out, 

I 

there  was  no  slipping  between  the  base  plate  and  the  support  blocks  at 
either  end  of  the  panel  and  so  the  back-up  block  was  not  used. 

d)  Loading  Head  and  Pivot 

The  analysis  assumed  that  the  panel  was  subjected  to  a  line  load 
applied  at  the  geometric  middle  surface.  Tf^re  were  two  difficulties  in 
achieving  this  ideal.  First,  it  was  physically  impossible  to  contact 
the  geometric  middle  surface.  A  line  load  acting  radially  from  above 
would  contact  the  outer  surface  of  the  panel.  This  would  result  in  a 
slight  offset  between  the  location  the  load  was  assumed  to  act  and  where 
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it  actually  did.  This  is  not  serious  until  the  cross  section  of  the 
panel  at  the  line  load  location  begins  to  rotate.  this  effect, 

coupled  with  the  load,  would  induce  a  small  moment  in  the  panel,  in 
addition  to  the  line  load.  This  could  ha \'e  been  accounted  for  in  the 
analysis  but  such  a  moment  was  deformat ion-^dependent,  making  the  govern¬ 
ing  equations  more  difficult  to  solve.  Thus  to  have  the  line  load  act 
at  the  geometric  middle  surface,  a  novel  loading  mechanism  was 
adapted.  Incorporated  into  this  mechanism  was  a  scheme  to  insure 
against  the  second  problem,  namely  that  no  moment  be  transferred  into 
the  panel  through  the  loading  fixture.  Figure  31  shows  the  loading  head 
and  pivot  fork  used  to  transmit  a  pure  line  load  to  the  geometric  middle 
surface  of  the  panel.  Figure  32  shows  the  pivoting  action  of  the  mecha¬ 
nism.  Three  small  bolts  through  the  panel  and  the  large  dowels  (shown 
in  fig.  24)  secured  the  arc-wise  location  of  the  line  load.  Centering 
the  panel  in  the  gap  between  the  large  dowels  guaranteed  that  the  line 
of  pivoting  passed  through  the  geometric  middle  surface.  Of  course  the 
pivots  themselves  prevented  any  moment  from  being  transferred  to  the 
panel.  Figure  33  shows  a  photograph  of  the  pivot. 

e)  Tongue  and  Connector 

The  pivot  fork  was  attached  to  the  bottom  end  of  two  long  steel 
plates.  The  long  plate  assembly,  or  tongue,  is  shown  in  fig.  34,  In 
the  analysis  it  was  assumed  the  line  load  moved  downward  radially  and 
deflected  the  panel,  with  the  simply  supported  ends  of  the  panel  remain¬ 
ing  stationary.  In  the  experiment,  the  vertical  displacement  of  the 


FIG.  51  Schematic  offloading  head  and  pivot 


panel  at  the  location  of  the  line  load  was  held  to  zero  and  the  simply 
supported  ends  moved  upward.  Due  to  the  assummed  shallowness  of  the 
panels  tested  this  was  essentially  the  same  as  holding  the  radial 
displacement  to  zero  at  this  point.  The  tongue  was  used  to  attach  the 
panel  to  the  stationary  upper  head  of  the  load  frame.  The  base  plate 
rested  on  the  lower  movable  head.  To  deform  the  panel,  the  lower  head 
was  moved  upward.  The  tongue  was  made  as  long  as  possible  for  the 
following  reason:  As  the  panel  deformed,  tangential  (u)  displacements 
resulted.  This  meant  that  the  line  load  position  moved  tangentially. 

In  the  analysis  it  was  assumed  that  the  line  load  acted  strictly  in  a 
radial  direction.  To  do  this  the  line  load  must  translate  tangentially, 
the  line  of  action  remaining  parallel  to  the  original  no-load 
direction.  Experimentally  it  was  very  difficult  to  translate  the  line 
load  tangentially.  An  infinitely  long  tongue,  fixed  at  the  top  and 
attached  to  the  panel  at  the  bottom,  would  satisfy  this  requirement. 
Practically  speaking,  a  long  tongue  is  all  that  was  possible.  The 
maximum  tongue  length  was  dictated  by  the  physical  dimensions  of  the 
load  frame,  A  connector  simply  attached  the  upper  end  of  the  tongue  to 
a  universal  joint  and  ultimately  the  fixed  upper  head.  Figure  35  shows 
another  overall  view  of  the  test  fixture  with  a^  panel  in  place. 

Instrumentation 

After  fabrication  the  specimens  Were  instrumented  with  350  ohm 
strain  gages  with  a  gage  length  of  0.25  in.  The  excitation  voltage  for 
the  strain  gages  was  2  volts.  The  strain  gages  were  mounted  back-to- 
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back  on  tne  specimens  so  that  the  middle  surface  strain  could  be  moni¬ 
tored.  The  location  of  the  gages  was  determined  by  the  amount  Of  load 
offset  for  the  specimen.  For  load  offsets  of  1S%  of  the  semi-opening 
angle  and  less,  the  gages  were  placed  at  the  quarter  points.  When  the 
load  offset  exceeded  15%  the  gages  were  placed  at  midspan  and  at  the 
quarter  point  furthest  from  the  load.  To  measure  the  radial  displace¬ 
ment,  dial  gages  mounted  radially  were  employed.  The  dial  gages  were 
capable  of  measuring  displacements  of  0.001  in.  and  greater.  Again  the 
location  of  the  dial  gages  depended  on  the  degree  of  the  load  offset. 
For  offsets  of  15%  and  less,  th^  dial  gages  were  located  at  the  quarter 
points  and  at  the  point  of  load  application.  For  offsets  greater  than 
15%  the  dial  gages  were  placed  at  the  point  of  load  application,  mid- 
span,  and  at  the  quarter  point  farthest  from  the  load.  The  tests  were 
carried  out  in  the  Engineering  Science  and  Mechanics  department  of 
Virginia  Tech  losing  a  displacement-controlled  Instron-loading  frame. 

The  load  was  determined  using  a  1000  pound  Instron  load  cell  placed  in 
siries  with  the  load  fixtiire.  The  straih  gage  data  was  acquired  using  a 
s)jst^m  consisting  of  Vishay  2^120  Wheatstone  bridges  and  signal  condi¬ 
tioner,  an  HP3495A  scanner  with  an  A/D  converter,  and  a  Tektronix  4051 
cotTiputer.  The  dial  gage  readings  were  manually  entered  into  the  Tektro¬ 
nix  computer  whi le  strain  gage  and  load  data  were  automatically  recorded 
on  digital  tape.  Figures  36  and  37  shows  the  experimental  set-up  with 
instrumentation  in  place. 
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FIG.  37  Photograph  of  experimental  set-up 

END  VIEW 
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Initial  Measurements 

Before  the  actual  testing  of  the  specimens  took  place,  some  initial 
measurements  related  to  panel  geometry  and  material  properties  were 
taken.  The  arc  length,  width,  and  thickness  were  measured  to  assure 
each  was  as  specified.  Also  the  positioning  of  the  holes  drilled  for 
load  positioning  were  checked  for  accuracyl  Tensile  tests  to  determine 
the  extensional  stiffness,  and  three  point  bend  tests  to  determine  the 
bending  stiffness  were  made.  The  tensile  tests  were  performed  on  strips 
cut  from  the  larger  panels  from  which  the  specimens  were  cut.  The  bend 
tests  were  conducted  on  the  actual  specimens.  Details  of  the  test 
procedure  and  results  are  contained  in  Appendix  C. 

The  initial  radius  of  each  specimen  was  measured  after  it  was 
fitted  into  the  test  fixture.  Measuring  the  radius  after  the  panel  was 
in  the  fixture  assured  that  the  measured  radius  was  indeed  the  initial 
radius  of  the  specimen  tested.  To  accomplish  the  measurement  of  the 

radilis,  a  dial  gage  and  a  traversing  mechanism,  capable  of  measuring  arc 

i  ^ 

length  location  to  within  0.05  in.,  were  used  to  obtain  the  spatial 

i  ' 

coordinates  of  several  points  on  the  mounted  specimen.  By  using  a 

1  , 

paf^abdlic  least-square  fit  to  the  data,  the  average  radius  of  the 

I 

specimen  was  determined.  Table  2  shows  the  measured  dimensions  for  the 
specimens  tested. 

Test  Procedure 

Having  made  the  initial  measurements  on  the  specimen,  excluding  the 
radius  measurement,  it  was  placed  loosely  in  the  test  fixture.  In  this 
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Table  2 


Measured  Test  Specimen  Geometry 


Specimen 

Number 

Radius 

(in.) 

Thickness 

(in.) 

1 

56.34 

.084 

2 

52.04 

.082 

3 

54.72 

.083 

4 

51.00 

.082 

5 

55.42 

.  .060 

6 

55.34 

.062 

50.35 

.063 

8 

50.09 

.063 

9 

52.56 

.042 

10 

55.41 

.043 

u 

57.33 

.042 

12 

64.7 

■  .082 

13 

62.4 

.081 

14 

63.8 

.081 

15 

50.41 

.022 

16 

— 

-- 

17 

-- 

-- 

18 

Note:  Of  specimens  15-18  only  specimen  15  was  actually  tested, 
thus  specimens  16-18  had  no  initial  measurements  made. 


position  the  strain  gages  vfere  connected  and  they  were  balanced  and 
zeroed.  The  si'ecimen  was  then  tightened  in  the  fixture  by  shding  the 
movable  support  block.  To  assure  a  tight  fit  the  specimen  was  slightly 
compressed  by  the  support  blocks  until  a  slight  bending  strain  was 
registered.  The  loading  head  was  then  attached  and  the  levelness  of  the 
fixture  and  the  vertical  orientation  of  the  tongue  were  checked  and 
corrected  as  needed.  The  Initial  radius  of  curvature  measurements  were 
then  taken.  After  this,  the  load  cell  reading  was  zeroed.  The  dis¬ 
placement  rate  was  set  at  0.01  in. /min.  and  readings  of  load,  strain, 
and  displacement  v/ere  recorded  at  every  0.01  In.  of  midspan  displace¬ 
ment.  To  record  the  dial  gage  readings,  the  loading  head  displacement 
was  stopped.  The  readings  were  then  entered  into  the  Tektronix.  During 
the  first  tests  the  support  blocks  were  monitored  with  dial  gages  for 
any  sign  of  elastic  or  rigid  body  displacement.  Since  the  blocks  showed 
no  sign  of  bending  or  slipping  relative  to  the  base  plate,  this  was 
discontinued  for  later  tests. 


Chapter  5 

CORRELATION  OF  EXPERIMENTAL  AND  THEORETICAL  RESULTS 

In  order  to  validate  an  analysis  it  is  necessary  to  compare  it  rfith 
experimental  results.  Here  the  analytic  results  presented  in  Chapter  3 
are  compared  to  experimental  results  obtained  using  the  apparatus  and 
test  fixture  described  in  Chapter  4.  So  as  to  make  the  comparison  with 
the  experiment  clearer,  the  analytic  results  shown  represent  only  those 
configurations  assumed  to  be  stable.  The  initial  measurements  made  on 
the  test  specimens  concerning  the  geometry  and  material  properties  were 
used  in  generating  the  analytical  curves  presented  with  experimental 
results.  Also,  the  small  initial  strain  in  the  panels,  due  to  fitting 
the  specimens  into  the  fixture,  has  been  subtracted  from  the  experi¬ 
mental  data  presented.  In  all  tests  performed  this  initial  strain  was 
less  than  50|ie  in  magnitude,  or  less  than  2,5X  of  the  maximum  strain. 
Both  the  experimental  and  analytic  results  are  presented  in  nondimen- 
sional  form,  with  the  exception  of  the  shape  data,  where  the  units  are 
inches. 

In  general,  the  agreement  between  theoretical  and  experimental 
results  was  good.  All  the  trends  predicted  by  the  analysis  were 
observed  in  the  experimental  results.  Although  the  agreement  was  seen 
to  be  good,  there  was  still  some  deviation  between  theory  and  experi¬ 
ment.  This  deviation  may  be  attributed  to  a  number  of  factors.  Perhaps 
the  most  significant  was  the  fact  that  the  specimens  were  found  to  have 
imperfect  shapes,  i.e.,  they  were  not  exactly  circular.  The  work  done 
by  Schreyer  [4]  shows  that  geometric  imperfections  have  much  the  same 
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effect  on  panel  .response  as  do  load  eccentricities.  Tnis  is  tne  most 
likely  reason  for  not  seeing  a  distinct  bifurcation  of  the  equilibrium 
path  When  the  load  was  exactly  at  midspan.  Other  possible  factors  for 
deviation  of  the  experimental  results  from  the  predictions  include  the 
possibility  of  friction-generated  moments  at  the  supports,  even  though 
great  care  was  taken  to  minimize  this.  The  friction  was  manifested  by 
small  load  jumps,  characteristic  of  end  support  sticking,  observed 
during  the  testing  of  the  panels.  Thickness  variations  observed  in  the 
test  specimens  during  the  initial  measurements  were  also  likely  to  cause 
differences  between  experiment  and  theory.  This  thickness  variation 
could  be  responsible  for  the  differences  between  the  predicted  and 
observed  maximum  surface  strains  seen  for  some  of  the  test  specimens. 

In  this  section  the  experimental  results  from  selected  test  cases 
are^resented.  These  cases,  along  with  their  associated  parameters,  are 
given  in  Table  3.  Since  the  trends  observed  for  most  of  the  panels 
tested  were  similar,  the  test  cases  shown  here  are  representative. 
Results  from  the  [90/0]s  panels,  which  were  quite  thin,  are  not  shown 
here.  Problems  with  these  specimens  developed  during  the  testing.  It 
appeared  that  friction  generated  at  the  supports  was,  at  times, 
sufficient  to  constrain  free  rotation  of  the  specimen.  The  specimen 
would  actually  bend  at  the  support  ends  and  so  the  specimen  behaved  at 
times  as  if  it  had  clamped  eHs  instead  of  simply  supported  ends. 

Figures  38-40  show  the  predicted  and  observed  response  for  a 
[{90/0)4]g  composite  specimen  with  a  12  in.  arc  length,  a  56.34  in. 
radius,  and  no  load  offset,  test  case  1.  The  analytic  predictions  are 
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Table  3 


Test  Cases 

Presented 

Test 

Case 

Number 

Specimen 

Number 

Lay-up 

Radius 

(in.) 

Arc 

Length 

(in.) 

Offset 

Angle 

{%) 

1 

1 

[(90/0)435 

56.34 

12 

0 

2 

2 

[(90/0)435 

52.04 

12 

15 

3 

3 

[(90/0)43s 

54.72 

12 

30 

4 

5 

[(9O/0333s 

55.42 

12 

.  0 

5 

6 

[(9O/0)33s 

55.34 

12 

15 

6 

9 

[90/45/0/-453s 

52.55 

12 

15 

7  . 

12 

[(90/0)43s 

64.7 

6 

0 

FIG.  33  Comparison  of  theory  and  experiment:  load-thrust 
^D  load-displacement  relations;  R=56.3d.  s=P  0 
r=0.n,  [(90/Oh].  " 


represented  by  the  solid  line  while  the  open  symbols  represent  the 
experimental  data.  The  analysis  predicts  that  such  a  panel  will  exhibit 
bifurcation  of  the  primary  equilibrium  branch.  The  load-displacement 
diagram  of  fig.  38  shows  that  the  observed  behavior  closely  follows  the 
predicted  response.  (The  darkened  symbols  in  fig.  38  correspond  to 
decreasing  values  of  thrust.)  From  the  discussion  of  bifurcation 
behavior  in  Chapter  3  this  suggests  that  the  initial  imperfections 
associated  with  the  panel  are  relatively  small.  The  load-thrust  diagram 
of  fig.  38  shows  more  clearly  that  the  actual  behavior  is  limit  point 
behavior.  Since  it  has  been  shown  that  geometric  imperfections  effect 
the  panel  response  in  much  the  same  wa^  as  do  load  eccentricities,  it 
would  be  beneficial  to  refer  to  the  discussion  in  Chapter  3  of  the 
effects  of  load  eccentricities  on  panels  exhibiting  bifurcation  behavior 
when  they  are  loaded  exactly  at  inidspan.  Appendix  1)  investigates  the 
effect  of  geometric  imperfections  on  the  panel  behavior.  It  should  be 
noted  here  that  the  inplane  stress  resultant,  N.,  is  assumed  to  be 
constant  (see  eq.  80).  Perhaps,  this  is  not  exactly  true  and  the  devia¬ 
tion  of  theory  and  experiment  is  due  to  a  slight  spatial  variation  of 
the  thrust.  To  experimentally  obtain  the  thrust  parameter,  q,  the  back- 
to-back  strain  gages  at  each  spatial  position  are  averaged  to  obtain  the 
middle  surface  strain.  The  thrust  parameter  is  then  calculated,  using 
eqs.  65,  105,  and  123,  at  each  point  and  is  then  averaged.  Despite  the 
symmetry  to  the  set-up,  for  a  given  specimen  and  load,  the  thrust  para¬ 
meter  at  each  strain  gage  location  did  not  have  the  same  numerical 
value.  This  suggests  that  q  may  not  have  been  spatially  constant. 
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Alternatively,  a  slight  misalignment  of  the  strain  gages  or  other  errors 
associated  with  strain  gage  measurements  could  have  made  a  very  accurate 
determination  of  q  impossible.  Also,  it  is  seen  from  the  load-disolace- 
ment  diagram  that  the  initial  panel  stiffness,  represented  by  the  slope 
of  the  equilibrium  path,  observed  experimentally  differed  from  that 
which  was  predicted.  This  difference  could  be  attributed  to  errors  made 
in  the  material  property  and  geometry  measurements  of  the  specimens. 

I 

The  parameter  X,  which  governs  the  panel  response,  is  sensitive  to  these 
measurements.  Appendix  E  explores  the  effect  of  the  initial  panel 
stiffness  on  the  overall  response.  Fig.  39  shows  the  panels  geometric 
configuration  of  various  displacement  levels.  The  panel  configurations 
predicted  by  the  analysis  are  purely  symmetric  up  to  the  bifurcation 
point.  After  this  the  configurations  become  asymmetric.  Though  not 
shown  in  fig.  39,  during  the  experiment  it  was  observed  that  the 
configurations  began  to  show  asymmetry  at  much  lower  load  levels  than 
were  predicted.  This  again  led  one  to  believe  that  geometric 
imperfections  were  present.  Figufe  40  shows  a  comparison  between 
theoretical  and  measured  strains.  The  maximum  surface  strains  predicted 
by  the  theory  are  somewhat  higher  than  the  observed  strains.  This  may 
have  been  due  to  local  thickness  variations,  since  only  an  average  panel 
thickness  was  used  in  the  analysis. 

Figures  41-43  show  the  response  for  a  [(90/0)435  12  in.  arcT  length 
composite  specimen  with  a  52.04  in.  radius  and  a  15%  load  offset,  test 
case  2.  With  no  load  offset  such  a  panel  shows  bifurcation  response. 
Ideally,  it  would  have  been  desirable  to  have  a  panel  with  a  radius  of 


I 


■( 
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FIG.  ^3  Comparison  of  theory  and  experiment:  surface^ 

STRAINS  vs.  applied  LOAD;  R=52.04^  S=12.0^  F=0.15, 


[(90/0)433 


curvature  equal  to  the  first  experimental  case  described.  Since  the 
initial  radius  of  curvature  was  quite  sensitive  to  the  span  between 
support  blocks,  it  was  difficult  to  obtain  the  same  radius  for  a  group 
of  specimens.  However,  since  radii  of  curvature  for  the  first  and 
second  experimental  cases  were  quite  close,  and  since  the  two  radii  of 
curvature  would  produce  bifurcation  behavior  in  the  no-offset  case,  it 
is  fair  to  say  that  the  only  difference  between  the  first  and  second 
experimental  cases  was  the  15%  load  offset.  This  second  experimental 
case  performed,  then,  as  expected.  A  comparison  of  this  offset  load 
case  and  the  previous  midspan  load  case  shows  a  decrease  in  the  experi¬ 
mentally  determined  limit  loads,  thus  verifying  the  prediction  that  such 
a  structure  is  sensitive  to  load  eccentricity.  The  load-thrust  and 
load-displacement  diagrams  of  fig.  41  show  good  agreement  between  thepry 
and  experiment.  The  same  is  true  of  the  shape  predictions.  The  :  urface 
strain  predictions,  shown  in  fig.  43,  deviate  from  the  observed  strains 
significantly. 

Figures  44-46  show  the  results  for  a  [(90/0)jj^]g  panel  with  a  3o% 

!  / 

offset.  The  panel  arc  length  Was  12  in.  and  the  radius  of  curvature  was 
54.72  in.  The  comparison  between  theory  and  experiment  for  the  load- 
thrust  relation  was  not  as  good  as  for  the  other  two  first  cases  but 
still  the  basic  trends  were  there.  For  offset  loadings,  the  maximum 
thrust  tends  to  increase  with  increasing  load  offset  for  panels  which 
exhibit  bifurcation  beiiavior  when  the  load  is  exactly  at  inidspan.  Thus 
the  force  on  the  supports  increases.  During  the  testing  of  this  speci¬ 
men  it  was  observed  that  friction  at  the  supports  had  more  of  an  effect 
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FIG.  46  Comparison  of  theory  and  experiment:  surface 

STRAINS  vs.  APPLIED  LOADj  R=54.72,  S=12.0>r’  =Q.3j 
[(90/0)^}^ 
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twr,  for  we  previous  tests.  The  frtctloo  was  detected  fro«  sltoht  load 
jumps  due  to  intermittent  sticking  of  the  supports,  thus  constraining 
rotation.  The  configurations  predicted,  fig.  45.  snoK  fair  agreement 
uitn  the  shapes  onserved.  The  predicted  and  observed  strains  once  again 

deviate  from  one  anothePc 

Figures  47-49  snow  the  response  for  a  [(90/0)33,  12  in.  arc  length 
panel  with  a  radius  of  55.42  in.  and  no  load  offset  while  figs.  50-52 
snow  the  response  of  a  very  similar  panel  (55.34  in.  radius)  with  a  15.. 
load  Offset.  For  the  midspan  loading  the  correlation  between  theory  and 
experiment  is  good,  particularly  for  the  load-displacement  relation  and 
the  strains.  Again  there  appeared  to  be  some  geometric  imperfections 
present  since  limit  point  behavior,  rather  than  bifurcation  behavior, 
was  observed  tor  the  specimen  with  no  load  offset.  This  was  also  born 

odt  by  the  observation  of  the  equilibrium  configurations  during  test¬ 
ing.  The  Shapes  observed  were  asymmetric  in  naUre  at  much  lower  load 
levels  than  were  predicted.  In  general .  however,  shape  correlation  was 
acceptable,  as  can  be  seen  in  fig.  4b.  The  offset  loading  case  also 
snows  extraordinary  agreement  between  theory  and  experiment.  Load- 
deflection.  load-thrust.  Shape,  and  the  surface  strain  predictions 

correlate  well  with  the  observed  results. 

The  results  for  a  [90/45/07-453,  12  f".  arc  length  panel,  with  a 
radius  of  52.56  in.  and  a  X5S  load  offset,  test  case  6.  are  shown  in 
fig.  53-55.  In  order  to  theoretically  predict  the  response  of  a  quasi 
isotropic  composite,  it  was  assumed  that  the  response  was  one-dimen¬ 
sional  and  that  the  only  nonzero  curvature  change  was  in  the  e  direc 
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FIG.  52  Comparison  OF  theory  and  experiment:  surface 
STRAINS  vs.  APPLIED  LOAD;  R=55.3^.  S=12.0. 

r=0.15.  [(90/0)J. 
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FIG.  5^  Comparison  of  theory  and 
CONFIGURATIONS;  R=52.56. 
t90/^5/0/-i}5]. 
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FIG.  55  Comparison  of  theory  and  experiment:  surface 
DRAINS  vs.  applied  LOADj  R=52.56^  S=12.0^ 
r=0.15.  [90A5/0/~/45]5 


117 


/ 


•v 


\z/ 


tion.  Tnus  a  raiuced  bending  stiffness  «s  used.  In  ream,.  Dj^  and 

025  Oeing  nonzero  Invalidate  tnis  assnnptlon.  Do,  to  the  constraint 
placed  on  the  test  specimen  by  the  fixture  it  «as  thought  that  these 
one-dimensional  assumptions  were  a  good  first  approximation.  A  compari¬ 
son  of  tneor,  and  experiment  snows  good  agreement.  The  correlation 
between  the  observed  and  predicted  load-displacement,  deformed  shapes, 
and  strains  are  also  seen  to  be  quite  good. 

Figures  56-57  show  the  response  of  a  [(SO/O),]^  6  in.  arc  length 
panel  with  a  radius  of  64.7  and  no  load  offset.  This  is  case  7  and  is 
the  last  case  discussed  here.  The  predicted  response  is  limit  point 
behavior  but  the  observed  response  is  that  of  monotonic  behavior.  The 
reason  for  this  discrepancy  between  theory  and  experiment  is  not 
Immediately  obvious.  Flexibility  of  the  supports  is  not  an  issue  here 
Since  the  forces  exerted  on  the  support  blocks  are  less  than  in  previous 
cases.  Error  in  the  initial  radius  measurement  has  also  been  ruled  out 
Since  it  was  double-checked  and  found  to  be  correct.  Also,  the  initial 
matfrtal  property  measurements,  recalling  their  effect  on  the 
parameter  x.  are  felt  to  be  accurate  enough  as  to  not  be  a  factor 
here.  It  seems  likely  that  there  is  a  "start-up"  problem  with  the 
experimental  set-up.  By  "start-up"  is  meant  low  load  and  tnr„,t  levels. 
I.e..  initial  portion  of  test.  This  could  be  the  result  of  mechanical  ’ 
Slack  or  dead  zones  in  the  fixture.  Becall  that  for  the  previous  test 
cases  there  was  a  discrepancy  between  the  predicted  and  observed  initial 
stiffnesses.  Appendix  E  ruled  out  the  possibility  of  error  in  initial 
measurements  as  the  cause  of  this  discrepancy.  Since  the  initial  mid- 
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FIG.  57  Comparison  OF  theory  and  experiment:  surface 

STRAINS  VS.  APPLIED  LOAD;  R=64.7^  S=6.0>r=0 

[(90/0)433 
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span  rise  of  the  panel  of  test  case  7  is  significantly  smaller  than  for 
the  other  test  cases,  it  is  possible  that  the  start  up  problem  effects  a 
much  greater  portion  of  panel  response,  thus  creating  the  poor  correla¬ 
tion  between  prediction  and  observation.  Whatever  the  cause,  it  was 
observed  in  tests  of  the  same  panel  with  load  offsets  of  15%  and  30%. 


Thus  the  phenomenon  was  repeatable. 

In  order  to  address  the  question  of  load  eccentricity,  the  limit 

point  loads  for  the  various  load  offsets  for  a  given  group  of  specimens, 
i.e.,  same  lay-up  and  arc  length,  were  plotted  along  with  an  upper  and 
lower  bound  on  the  theoretical  predictions.  Here  the  upper  bound  corre¬ 
sponds  to  the  curve  for  the  panel  with  the  smallest  value  of  while 
the  lower  bound  corresponds  to  the  curve  for  the  panel  with  the  largest 
value  of  X.  Since  the  radius  of  each  panel  within  the  group  changes,  it 


is  necessary  to  plot  the  upper  and  lower  bound.  Figure  58  shows  the 
(Jlot  of  limit  load  vs.  percent  load  offset  for  the  [(90/0)435  i2  in.  arc 
ength  specimens.  From  this  figure  it  seems -that  the  sensitivity  to 
oad  eccentricities  predicted  is  indeed  close  to  that  which  was 
(Observed.  Similarly  fig.  59  shows  the  sensitivity  plot  for  the 


[(90/0)335  12  in.  arc  length  specimens,  and  fig,  60  shows  the  sensi¬ 
tivity  for  the  [90/45/0/-4535  specimens.  It  appears  the  observed  sensi 
tivity  has  been  accurately  predictedo 


FIG.  58  Comparison  of  theory  and  experiment:  sensitivity 
OF  buckling  load  to  load  eccentricity^  [(90/0)^] 

SPECIMEN  GROUP 
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FIG.  59  Comparison  of  theory  and  experiment: 

OF  BUCKLING  LOAD  TO  LOAD  ECCENTRICITY^  [(90/0)31 
SPECIMEN  GROUP 


Fi6.  60  Comparison  of  theory  and  experiment:  sensitivity 

OF  BUCKLING  LOAD  TO  LOAD  ECCENTRICITY^ 
[90A5/0/-45U  SPECIMEN  GROUP 


Chapter  6 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  study  presented  In  this  thesis  was  intended  to  extend  the 
knowledge  of  the  static  behaviorof  shallow  orthotropic  curved  panels. 
The  principle  aims  of  the  study  were  to  determine  the  effects  of  ortho¬ 
tropic  material  properties  on  the  behavior  of  shallow  curved  panels,  to 
experimentally  verify  the  analysis,  and  thus  ultimately  lay  the  ground 
work  for  a  dynamic  analysis  of  the  same  problem.  An  exact  solution, 
taking  into  account  load  eccentricity,  was  developed  for  a  simply 
supported  circularly  cylindrical  curved  panel  subjected  to  a  line  load 
along  a  cylinder  generator.  Also,  an  experimental  apparatus  was 
designed  to  load  the  curved  panels.  The  analytical  predictions  were 
compared  with  experimental  measurements. 

From  this  study  the  following  conclusions  may  be  made: 

1.  Analytically,  three  types  of  panel  behavior  were  seen  for 
panels  loaded  exactly  at  midspan;  monotonic,  limit  point,  and 
bifurcation.  Only  two  types  of  behavior  were  exhibited  when 
the  load  was  offset;  monotonic  and  limit  point.  The  type  of 

behavior  predicted  depends  solely  on'  one  parameter,  X.  Geome- 

\ 

try,  material  properties,  and  stacking  sequence  influence  X  for 
an  orthotropic  panel.  For  an  isotropic  panel  X  is  affected  by 
geometry  only. 

2,  Those  panels  which  exhibit  bifurcation  behavior  when  loaded 
exactly  at  midspan  were  seen  to  be  imperfection-sensitive. 

This  was  demonstrated  both  analytically  and  experimentally. 


InttUl  geometric  Imperfections  were  detected  end  guantifled 
«nile  making  initial  radius  measurements.  Due  to  tne  imperfec¬ 
tions  present,  innerent  of  any  real  structure,  bifurcation 

benavior  was  not  seen.  However  tne  observed  response  did 

Closely  approximate  bifurcation  behavior. 

The  overall  correlation  between  theory  and  experiment  was 
good.  This  implies  that  tne  solution  metbod  used  is  reasonably 
accurate.  Predictions  of  tne  limit  loads  were  guite  accurate. 
The  predictions  of  strain  in  tne  panels  were  fair  for  tne  most 
part.  However,  a  spanwise  spatial  variation  in  the  middle 
surface  strain  was  observed  during  testing.  This  is  counter  to 
tne  assumption  tnat  tne  strain  was  constant  along  tne  span. 

Due  to  tne  small  bending  stiffness  innerent  to  thin  laminates, 
it  was  impossible  to  obtain  accurate  results  for  these 
laminates  with  tne  current  experimental  set-up.  The  friction 
related  moments  generated  at  tne  ends  of  tne  simply  supported 
panels  at  times  caused  these  panels  to  behave  as  1f  they  were 

clamped. 

.  in  analyzing  tne  test  results  there  seetas  to  nave  been  "start- 
uj)“  or  initial  load  and  deformation  problems  with  the  experi¬ 
mental  set-up.  This  is  borne  out  by  the  lack  of  correlation 
between  prediction  and  observation  for  the  low  thrust 
regimes.  This  would  explain  the  discrepancies  in  initial 
stiffnesses  as  wel’.  as  the  total  lack  of  correlation  for  test 
case  7  in  which  the  maximum  thrust  was  expected  to  be  rela¬ 


tively  small. 
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As  a  result  of  this  work,  a  number  of  topics  can  be  singled  out  for 
further  investigations.  Of  these  topics,  a  number  are  relatively  minor 

extensions  of  the  previous  work.  These  are: 

1.  Modify  the  analysis  to  include  different  boundary  conditions; 

e.g.,  ffictional  effects. 

2.  Include  the  effects  of  geometric  imperfections  in  the  analysis. 

3.  Conduct  further  static  tests  with  various  material  properties 
and  geometries.  By  varying  the  material  properties  the  effects 
of  orthotropicity  described  in  this  text  may  be  experimentally 
verified. 

4.  Investigate  the  “start-up"  problem  encountered  during  testing. 

In  addition  to  the  above,  recommendations  requiring  significant 
extensions  of  this  work  are: 

1.  Extend  the  analysis  to  include  the  effects  of  shear  deformation 
and  to  take  into  account  angle  ply  laminates.  This  would  allov; 
a  far  greater  number  of  laminates  to  be  investigated  and  would 
perhaps  be  an  intermediate  step  to  the  dynamic  analysis. 

2.  Extend  the  analysis  into  the  dynamic  regime.  Ultimately  it  is 
hoped  that  an  understanding  of  the  behavior  of  a  composite 
fuselage,  upon  impacting  a  flat,  rigid  surface,  may  be 
obtained.  By  doing  so  the  crashworthiness  of  a  composite 
fuselage  may  be  assessed. 


V 
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Appendix  A 

COEFFICIENTS  OF  LOAD-THRUST  RELATIONSHIPS 
The  solution  technique  applied  to  the  governing  differential  equa¬ 
tions  discussed  in  Chapter  3  requires  the  solution  of  a  quadratic  equa¬ 
tion.  This  equation  relates  either  the  applied  load  or  a  constant  of 
the  general  solution  of  ^{V)  to  the  thrust.  Here  the  functional  forms 
of  the  coefficients  to  the  quadratic  equation  are  presented.  For  each 
case  the  coefficients,  found  to  be  transcendental  functions  of  the 
thrust,  are  found  by  substituting  the  functional  form  of  the  radial 
displacement  Into  the  integrated  strain-displacement  relationship, 
expressed  nondimensionally  by  eq,  125. 


Ml dspan  Loading:  q  /  0  and  q  /  fi  =  1,  3,  5,  ... 


The  load-thrust  relationship  is  of  the  form. 


Ac'^  +  B5  +  C  =  0  , 


when  q  nn,  n  =  1,  2,  3,  ...  •  Here  A,  B,  and  C  are, 


^  ^ _ 3_  ^  tan^(q)  3tan(q 

16q^  16q‘^  16q^ 


B  =  1  .  -tanM- 


r  -  _  i  +  1  -  STn(.2_^L- 

'r?  ®  4q^COS^(q)  8q^COS^(q) 
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When  q  =  nu,  n  =  1,  2,  3,  ...  the  integrated  strain-displacement  rela¬ 
tionship  yields. 


+  C  =  Db^. 


This  relationship  relates  the  applied  load  to  the  magnitude  of  the 
asymmetric  displacement.  The  coefficients  are, 


B  [1  -  (.1)"] 

c  -  1  jai  1 


4q 


Midspan  Loading: 


D=a:. 

4 


n  —  1,  3,  5,  ... 


The  integrated  strain-displacement  relationship  yields. 


a(aj2  +  BA%  C  =  0  , 


to  determine  a  constant,  A„,  from  the  general  solution  of  6  (r)  given  by 
eq.  134.  Here  the  coefficients  are. 


(All) 
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(A12) 


131/ 


?■*  4,^  q- 


(AU) 


Fr.r.entric  Loading:  q  0  and  q  n  =  1.  2^ 


The  load-thrust  relationship,  when  load  eccentricity  is  present  is. 


+  C  =  0 


(A14) 


For  the  offset  load  case  the  coefficients  of  the  quadratic  equation  are 
functions  of  the  load  offset  as  well  as  the  thrust.  The  coefficients 

are: 


9  — 

3  (r  "r  sin(2qF  )  _  sin  (2qr  ) 
^  "  leq'^  8q^  16q^sin(2q)  16q^tan(2q) 


sin^(qr  ) 


cos^lqr  )sin(2q)  ^  _ 

^  I6q^  4q^tan(2q)sin(2q) 


(A15) 


3cos^(qf  )tan(q)  _  sin^(2qr  )tan(ql  _  s1n^(qrj. 
16q^  16q°  ^  tan(q) 

+  5sin^(qr) 

16q^tan(q) 

„  r*sin(qr)  tanqcos(qr)  _  ^  _cosj;qrj_ 


4q^cos(q) 


4q^C03(q) 


2q^cos (q)  2q^ 


1 


c  +  — o — ^ 

6  8q3cos'^(q)  4q'^cos^(q) 


{A16) 

(A17) 
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Eccentric  Loading:  q  =  ^»  ^  3» 

For  both  odd  and  even  values  of  n,  the  integrated  strain 
ment  relationship  yields, 

A(A*)^  +  8A*+  C  =  0. 

'  n  n 


For  odd  values  of  n  the  coefficients  are: 


B  =£££iLa?:J.  (.1)"  (f  .1)  sin^(qr  ) 

q  q 


+  Iq  (-1)"  [cos(2qr  )  -  cos(2q)] 


"  4q^  2qcos{qr  ) 

.^.i!£!|3Li]  .[(7)2-  1)  (-1)" 

2q^  2q^ 

[cos(qr  )  +  1]  +  ii-.). A''”  (y  )  (-1)'^ 
2q'^cos'^(qr  ) 

_  .si^^(qr  2  ,  cos(qr  ) 

2q^C0S^(qr  )  q^ 


r  sin(qr  )  sin^fqr  ) 


q  cos(qr  ) 


displace- 


(A18) 


(A19) 


(A20) 


(A21) 


i 
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For  even  values  of  n  the  coefficients  are: 


{A22) 


(A23) 


c  =i--a- 


(A24) 


i 

;  ,  •  1  ■ 

.{■  ’  ‘  '  i  ' 

/. 

»• 

\TJ 

'-Auv; 

Of'**'  '***•“  4  4*X/ 

OF  FOOi-1  CUAU3“V 

Appendix  B 

EFFECT  OF  ORTHOTROPY  ON  PANEL  BEHAVIOR 

The  parameter  X  is  defined  as. 

/ 

4  _  A22R‘‘b'’'  . 

A.  . .  » 

(Bl) 

^22 

where  A22»  the  extensional  stiffness,  and  0224  the  bending 

stiffness, 

are  given  by. 

- 

* 

^22  ""  ^22^k^^k"  ^k-1^ 

(B2) 

^22  I  ^^22^k^^k"  ^k-1^  * 

(B3) 

Concentrating  on  cross-ply  laminates,  where  each  lamina  is 

composed  of 

the  same  material,  the  reduced  transformed  stiffness  is. 

Q22,  ‘i^or  0°  laminae 

^22  " 

(B4) 

Ql|,  for  90°  laminae. 

The  Q's  are  the  reduced  stiffness  and  are  defined  as. 

/ 

/ 

• 

Ep9 

^22  ~ 

^-'’12''21 

(B5) 

Q  -  ..^11 

Wii  “  > 

1-V12V21 

(B6) 

. 
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wnere  is  the  Young's  modulus  in  the  fiber  direction  and  E22 
Young's  modulus  perpendicular  to  the  fiber  direction.  By  definition  it 
i  s  assumed  that  E^,  ^  >  E22. 

Using  the  above  definitions  for  the  stiffnesses,  a  general  expres¬ 
sion  for  the  extensional  stiffness  and  the  bending  stiffness  can  be 
developed  in  terms  of  the  number  of  laminae,  the  stacking  sequences,  and 
the  elastic  moduli  of  the  laminae.  For  convenience  these  expressions 
are  shown  only  for  the  [(gO/O)^]^  family  of  laminates.  For  these  lami¬ 
nates,  the  extensional  and  bending  stiffnesses  are. 


j-  (E^j+  E22) 

^22 "  1 

I-V12V21 

2  3 

[(K+3)Ejj+  (N-3)E22] 

O22 "  ;  * 

1  -  V12V21 


(87) 


(B8) 


where  I'l  is  the  number  of  laminae  (N=4n)  and  t  is  the  thickndss  of  each 
lamina.  Substituting  eqs.  B7  and  B8  into  eq.  B1  yields. 


where  use  is  made  of  the  fact  S  =  2Rp. 

To  emphasize  the  importance  of  taking  the  orthotropy  of  the  mater¬ 
ial  into  account,  three  shallow  circularly  cylindrical  panels  are  com- 
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pared  as  to  the  behavior  exhibited.  Each  panel  is  assumed  to  have  the 
same  geometry,  only  the  material  used  In  the  fabrication  of  the  panels 
is  varied.  For  the  comparison  aluminum,  glass-epoxy,  and  graphite-epoxy 
panels  are  used.  Table  B1  summarizes  the  geometry,  material  properties, 
and  calculated  values  of  X  for  each  of  the  panels.  The  value  of  \  for 
the  aluminum  panel  is  such  that  bifurcation  on  the  primary  equilibrium 
path  would  occur  when  the  panel  is  loaded  exactly  at  midspan.  However, 
when  glass-epoxy  or  graphite-epoxy  is  used,  limit  point  behavior  would 
be  expected.  Thus  by  increasing  the  degree  of  orthotropy,  the  panel 
behavior  changes  from  bifurcation,  where  a  asymetric  response  would  be 
expected  when  in  displacement  control,  to  limit  point  behavior,  where  a 
jump  in  load  and  thrust  would  occur. 

For  a  composite  material,  the  number  of  laminae,  as  well  as  the 
degree  of  orthotropy,  effect  the  value  of  X.  Figure  B1  shows  the  effect 
of  the  number  of  laminae,  as  well  as  the  degree  of  orthothopy  present, 
on  the  calculated  value  of  X.  With  an  increasing  number  of  laminae,  the 
effect  of  orthotropy  becomes  less  pronounced. 

Stacking  arrangement,  though  not  explored  specifically  here,  could 
be  important  in  determining  X  and  thus  panel  response.  Though  they 
would  have  the  same  weight,  two  panels  each  with  16  laminae  could  be 
made  to  behave  quite  differently,  depending  on  how  the  laminae  are 
arranged. 
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Table  B1 

Panel  Geometries  and  Material  Properties 


Aluminum 

Panel 

Glass-Epoxy 

Graphite-Epoxy 

Arc  Length  (in.) 

6 

6 

6 

Radius  (in.) 

45 

45 

45 

Thickness  (in.) 

0.08 

0.08 

0.08 

Lay  up 

— 

[904/0435 

1:904/0433 

Ell  (PS’) 

10  X  10® 

7.8  x  10® 

25  X  10® 

E22  (psi ) 

10  X  10® 

2.6  X  10® 

1  X  10® 

E11/E22 

1.0 

3.0 

25.0 

X 

2.94 

2.72 

2.58 
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AND  LAY-UP  ON  PANEL  RESPONSE 


Appendix  C 

determination  of  EXTENSIONAL  STIFFNESS  AND 
BENDING  STIFFNESS 

As  mentioned  in  Chapter  4.  the  initial  measurements  made  on  the 
test  specimens  included  the  determination  of  the  extensional  and  bending 
stiffnesses.  To  accurately  compare  the  results  of  the  analysis  and  the 
experiment  it  is  necessary  to  have  correct  material  properties  to  input 
into  the  analysis.  Here  the  test  set-ups  and  procedures  used  for  the 
measurement  of  both  the  extensional  stiffness  and  bending  stiffness  are 
described. 

Extensional  Stiffness 

The  extensional  stiffnesses  for  each  layup  were  obtained  from 
tensile  test  data.  A  tensile  test  specimen  was  machined  from  each  of 
the  large  panels,  as  were  the  curved  panel  test  specimens.  The  tensile 
coupon,  0.5  X  12  in.,  was  cut  from  the  panel  such  that  the  stiffness  in 
thd  X  direction  was  actually  measured.  Since  the  specimens  were  cross- 
ply  laminates,  the  stiffnesses  in  the  x  direction  and  0  direction  should 
hive  been  equal.  The  specimens  v^ere  tested  in  the  Instron  Tensile 
Machine  using  an  Instron  load  cell  to  measure  applied  load.  Back-to- 
back  strain  gages  were  mounted  at  the  center  of  the  gage  section.  The 
strain  gages  used  were  120  ohm  gages  with  approximately  4  feet  of  lead 
wire.  The  excitation  voltage  for  the  strain  gages  was  2  volts.  The 
strain  gage  and  load  data  were  acquired  by  using  a  digital  data  acquisi¬ 
tion  system.  For  each  test  performed  the  slope  of  the  load  vs.  strain 
curve  was  calculated  using  a  least  squares  curve  fit.  This  slope. 
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proportional  to  the  laminate  Young's  modulus  in  the  x  direction,  was 
then  divided  by  the  width  of  the  tensile  specimen  to  obtain  the 
extensional  stiffness.  Each  test  was  performed  twice  in  order  to  obtain 
results  as  accurately  as  possible.  The  results  for  the  layups  investi¬ 
gated  are  summarized  in  Table  Cl. 

Bending  Stiffness 

The  bending  stiffnesses  for  each  of  the  curved  panel  specimens 
tested  were  obtained  by  performing  three-point  bend  tests  on  the  speci¬ 
mens  themselves.  The  specimens  were  placed  on  the  fixture  so  that  the 
panel  was  simply  supported,  restricting  end  displacement  in  the  vertical 
direction  while  allowing  freedom  to  displace  horizontally.  The  load  was 
applied  through  a  hanger  located  at  midspan  to  which  calibrated  dead 
weights  were  attached.  A  schematic  of  the  test  fixture  used  for  the 
three-point  bend  tests  is  shown  in  fig.  Cl.  A  dial  gage  Capable  of 
resolving  displacements  to  0.001  in.  was  used  to  measure  the  vertical 
displacement  at  midspan.  Load  vs.  displacement  data  was  recorded  and 
the  slope  determined  using  a  least-squares  fit  to  the  data.  Assuming  a 
linear  straight  beam  theory  was  suitable  for  describing  the  behavior  of 
the  bend  specimen,  the  bending  stiffness  was  obtained  from  the 
expression, 

D 

^22  48w 


(Cl) 


X 

^  J 

t  / 

. 

/ 

/' 

i  i! 

1  ^  1 

■' 

■  / 
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Table  Cl 

Measured 

Panel  Stiffnesses 

« 

• 

Test  No. 

Lay-up 

A22  (Ib/in) 

D22  (lb-in) 

i 

1 

[(90/0)4]s 

7.8194  X  10^ 

616.2 

i  i 

2 

[(90/0)4]s 

7.8194  X  105 

606.3 

] 

3 

C(90/0)4]s 

7.8194  X  10^ 

617.2 

i 

4 

[(90/0)4]s 

7.8194  X  10^ 

612.7 

;  f 

i 

5 

C(90/0)3^s 

5.85  X  10^ 

.  263.5 

;  !  ■ 

6 

[{90/0)3]s 

5.85  X  10^ 

265.4 

1  1 

> 

7 

[(90/0)3]s 

5.85  X  10^ 

268.2 

/ 

8 

[(90/0)3]s 

5.85  X  105 

260.1 

•  1 

9 

C90/45/0/-45]s 

3.175  X  10^ 

90.9 

} 

i  i  ‘ 

10 

C90/45/0/-45]s 

3.175  X  105 

90.5 

'  i 

:  \ 

11 

[90/45/C/-45]s 

3.175  X  10^ 

93.3 

■» 

12 

[(90/0)4]s 

7.8194  X  105 

613.5 
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where  p  is  the  load  intensityj  w  is  the  mi dspan. displacement  and  L  is 
the  span  between  supports.  The  results  of  the  three-point  bend  tests 
are  summarized  in  Table  Cl. 
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Appendix  D 

EFFECT  OF  GEOMETRIC  IMPERFECTIONS  ON  PANEL  BEHAVIOR 

As  was  mentioned  earlier,  Schreyer  [4]  found  the  effects  of  geo¬ 
metric  imperfections  on  the  behavior  of  isotropic  curved  panels  to  be 
qualitatively  the  same  as  the  effects  of  load  offsets.  The  geometric 
imperfections  investigated  in  Schreyer's  study  were  of  the  same  form  as 
the  asymmetric  deformation  mode  of  the  curved  panels.  In  the  present 
study  it  was  thought  that  similar  geometric  imperfections  were  partially 
responsible  for  some  of  the  deviation  between  theoretical  and  experi¬ 
mental  results  given  in  Chapter  5.  Since  the  location  of  the  holes 
drilled  into  the  panel  fcr  locating  the  loading  head  were  measured  and 
felt  to  be  accurately  centered,  load  eccentricity  v/as  not  felt  to  be  a 
problem. 

Figures  38-40  show  the  predicted  and  observed  response  for  a 
C(90/0)4]s  composite  specimen  with  no  load  offset.  The  behavior  expec¬ 
ted  of  this  panel  was  bifurcation  behavior.  While  the  load-displacement 
diagram  of  fig.  38  appears  to  show  bifurcation  behavior,  the  load-thrust 
diagram  leads  one  to  believe  that  some  imperfection  may  be  present.  For 
this  test  case  a  number  of  measurements  were  made  along  the  entire  span 
of  the  specimen  for  determining  the  average  radius.  (For  subsequent 
tests  these  measurements  were  made  only  over  half  the  span.)  From  these 
measurements  it  was  possible  to  determine  the  deviation  of  the  initial 
configuration  from  the  assumed  perfect  cylindrical  configuration. 

Figure  Dl  shows  the  measured  deviation  of  the  initial  configuration 
from  a  perfect  cylindrical  geometry.  The  perfect  configuration,  yo(x). 
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is  assumed  to  be  a  cylinder  with  a  radius  of  curvature,  which  is 
determined  by  a  least-squares  fit  of  the  data  taken  along  the  span  of 
the  curved  panel.  The  measured  configuration  is  denoted  by  y(x).  From 
this  figure  it  is  seen  that  the  deviation,  or  imperfection,  tends  to 
taka  on  the  form  of  the  asymmetric  deformation  mode  discussed  earlier. 
Since  geometric  imperfections  were  found  from  the  initial  measurements 
and  were  found  to  be  of  the  same  form  as  the  asymmetric  deformation 
mode,  as  was  the  case  in  Schreyer's  work,  it  was  a  distil.  '■  possibility 
that  some  of  the  differences  betv/een  theory  and  experiment  are  due  to 
these  imperfections. 

To  explore  this  point  further  a  5%  load  offset  was  compared  with 
the  theoretical  midspan  loading  and  the  experimental  points.  The  5% 
load  offset  was  used  to  show  the  qualitative  effect  of  a  small  imperfec- 
,tion  in  the  initial  configuration.  Figure  D2  shows  the  theoretical 
curves,  both  a  0%  and  a  5%  load  offset,  and  experimental  data  for  the 
load-thrust  and  load-displacement  relations  for  the  panel  of  test  case  1 
(see  Chapter  5).  The  5%  offset  case  tends  to  more  closely  approximate 
the  data  taken  from  the  test  over  a  greater  portion  of  the  response 
curves. 

Finally,  fig.  D3  shows  the  spatial  shapes  for  the  first  18  load  or 
displacement  points  of  fig,  02,  The  predicted  shapes  are  of  course 
syrrenetric  with  respect  to  midspan.  For  the  first  6  or  7  data  points, 
the  measured  shapes  are  also  symmetric.  However,  the  8th  point  shows 
that  the  actual  shape  has  lost  symmetry,  the  left  side  deflecting  more 
than  the  right  side.  At  points  13  and  14  the  asymmetry  has  shifted,  the 


FIG.  D2  Comparison  of  theory  and  experiment; 

LOAD-THRUST  AND  LOAD-DISPLAC^ENT 

relations;  R=56.3^.,  S=12.0^  P=0.0  and 
r=0.05.  C(90/0)J  ■ 


FIG.  D3  Comparison  of  theory  and  EXPERiMEjiT:  Epuilibrium 
configurations;  R=56.34^  S=12.0^  F  =0.0^ 

[(90,0) 
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right  side  deflecting  more.  These  deviations  in  synmetry  of  the  spatial 
shape  are  felt  to  be  due  directly  to  initial  imperfection  of  fig.  01. 

Even  though  the  work  presented  in  this  appendix  is  only  for  one 
specific  case,  in  all  likelihood  some  form  of  geometric  imperfection  is 
present  in  each  of  the  other  test  cases  as  well.  The  results  shown  here 
show  that  the  geometric  imperfection  found  was  of  the  appropriate  form 
to  cause  an  imperfect  response  and  that  by  taking  the  initial  imperfect 
configuration  into  account,  the  predicted  response  more  closely  models 
the  panel  behavior. 


Appendix  E 

EFFECT  l)F  INITIAL  STIFFNESS  ON  PANEL  BEHAVIOR 
In  Chapter  5  it  was  noted  that  the  initial  stiffness,  i.e.,  the 
initial  slope  of  the  load-displacement  relation,  predicted  by  the  theory 
and  the  observed  initial  stiffness  did  not  correlate  well.  Perhaps  some 
of  the  discrepancy  between  the  theory  and  experiment  was  due  to  this 
initial  slope  difference.  Here  the  effect  of  the  initial  stiffness  on 
the  overall  response  is  explored  with  the  expectation  that,  by  correctly 
modeling  the  linear  behavior,  the  nonlinear  behavior  can  be  closely 
modeled.  The  initial  panel  stiffness  may  be  determined  by  examining  the 
linear  response  of  the  panel.  In  order  to  determine  this  stiffness,  tne 
theorems  of  Castigliano  are  used  for  the  no-load  offset  case.  A  full 
investigation  of  this  is  found  in  [11]. 


The  complimentary  strain  energy  for  the  panel  is. 
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R  de. 
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For  a  shallow  panel  the  statically  admissible  forces,  found  from  a  free 
body  diagram,  are 

Ng  =  -T  {E2) 

where  T  is  the  horizontal  support  reaction.  Substituting  these  forces 
into  eq.  El  yields 
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V  =  ^  [C /  -  c  TP  +  C  P*^] 
^22 


(E4) 


in  which 


^  *^22  .  2R^p^ 

c  =^5!§i 

•P  24 


=  5V 
12 


(E5) 

(E6) 

(E7) 


Applying  Castigliano's  second  theorem  to  the  complimentary  strain  energy 
determines  the  horizontal  reaction  in  terms  of  the  applied  load.  Then, 
by  applying  Castigliano's  first  theorem  to  the  complimentary  strain 
energy  an  expression  for  the  initial  stiffness  is  obtained.  Thf  initial 
St. ff nest  is. 


K  = 


'22 


2Ci 


Rp  4CyCp-Cjp 

Substituting  for  the  constants  Cy,  Cjp,  and  Cp  the  stiffness  is 


rewritten  as. 


K  = 


16D,, 

(-#■) 


1  375 

T  ■  m 

where  use  was  made  of  the  relationships. 


7  ^ 


(E8) 


(E9) 


S  =  2Rp 


(ElO) 
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■  Froni  eq.  E9  it  is  seen  that  the  geometric  and  material  properties 
of  the  panel  govern  the  initial  response.  In  order  to  model  the  non¬ 
linear  response  accurately  it  seems  that  the  predicted  linear  response 
should  match  the  observed  response.  This  may  be  accomplished  by  'jugg¬ 
ling'  the  panel  parameters  in  the  theoretical  formulation  until  the 
predicted  response  matches  that  of  the  observed.  Table  El  shows  the 
effect  of  changing  the  panel  parameters  on  the  initial  stiffness  of  the 
panel.  Each  column  shows  the  effect  of  a  percentage  decrease  of  the 
parameters.  The  initial  parameters  used  are  those  of  test  casd  1  from 

Chapter  5. 


Upon  examining  the  predicted  and  observed  initial  stiffnesses  it  is 
seen  that  the  parameters  must  be  changed  substantially  to  match  the 
stiffness.  By  changing  the  panel  parameters  ko  dramatically,  the  pre¬ 
dicted  character  of  the  nonlinear  response  is  changed  considerably. 

This  is  shown  in  fig.  El.  In  this  figure,  the  initial  stiffness  in  the 
analysis  is  adjusted,  through  \,  to  match  the  observed  initial 
stiffness.  It  is  seen  that  although  the  initial  portions  of  the  load- 
displacement  and  load-thrust  agree,  the  overall  comparison  is  very 
poor.  Thus  it  seems  inaccuracies  in  initial  stiffness  could  not  be 
responsible  for  the  differences  between  experiment  and  theory  observed 
in  this  study.  Also,  since  small  inaccuracies,  say  5%,  in  the  panels 
geometric  and  material  properties  have  such  a  small  impact  on  the 


t 

1 


V' 


: 

154 
Table  El 

Effect  of  Panel  Parameters  on  Initial  Stiffness 


K/Ko 


1% 

5% 

loj; 
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1.008 
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D22 
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.968 

,ps 

*22 
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.981 

.951 
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U015 
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1.152 
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initial  stiffness,  it  seems  unlikely  t.,dt  poor  measurements  of  toe 
panel's  geometric  and  material  properties  would  be  responsible  for 
differences  between  observed  and  predicted  panel  response. 


